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Abstract 

We examine in detail the relative equilibria in the four-vortex problem where 
two pairs of vortices have equal strength, that is, T\ = T2 = 1 and T% = = m 
where m E R — {0} is a parameter. One main result is that for m > 0, the convex 
configurations all contain a line of symmetry, forming a rhombus or an isosceles 
trapezoid. The rhombus solutions exist for all m but the isosceles trapezoid case 
exists only when m is positive. In fact, there exist asymmetric convex configurations 
when m < 0. In contrast to the Newtonian four-body problem with two equal pairs 
of masses, where the symmetry of all convex central configurations is unproven, the 
equations in the vortex case are easier to handle, allowing for a complete classification 
of all solutions. Precise counts on the number and type of solutions (equivalence 
classes) for different values of m, as well as a description of some of the bifurcations 
that occur, are provided. Our techniques involve a combination of analysis and 
modern and computational algebraic geometry. 
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1 Introduction 

The motion of point vortices in the plane is an old problem in fluid mechanics that was 
first given a Hamiltonian formulation by Kirchhoff [T7j. This widely used model provides 
finite-dimensional approximations to vorticity evolution in fluid dynamics. The goal is to 
track the motion of the point vortices rather than focus on their internal structure and 
deformation, a concept analogous to the use of "point masses" in celestial mechanics. As 
with the Newtonian ra-body problem, an important class of homographic solutions exist 
where the configuration is self-similar throughout the motion. Such solutions are described 
as stationary by O'Neil [2Tj and are not limited to just relative equilibria (rigid rotations), 
but also include equilibria, rigidly translating solutions and collapse configurations. 

In this paper we focus primarily on the relative equilibria of the four-vortex problem 
when two pairs of vortices have the same vortex strength. Specifically, if ^6 1- {0} 
denotes the vortex strength of the z-th vortex, we set i\ = T 2 = 1 and T 3 = T 4 = m, 
treating m as a real parameter. Our main goal is to classify and describe all solutions as m 
varies. Four-vortex configurations which are not collinear (nor contain any three vortices 
which are collinear) can be classified as either concave or convex. A concave configuration 
has one vortex which is located strictly inside the convex hull of the other three, whereas 
a convex configuration does not have a vortex contained in the convex hull of the other 
three vortices. 

The symmetry and asymmetry of solutions plays a major role in our analysis. In fact, 
part of the motivation behind our work was to determine whether symmetry could be 
proven for this choice of vorticities when m > and the configuration is assumed to be 
convex. This question, while solved in the Newtonian four-body problem when the equal 
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Shape 


m G (-1, 11 

V 'J 


Type of solution (number of) 


Convex 


m — 1 


Square (6) 




< m < 1 


Rhombus (2), Isosceles irapezoid (4) 




— 1 < m < 


Rhombus (4) 
Asymmetric (8) 




-1/2 < m < 


Kite^a (4) 




m* < m < — 1/2 


Kite 12 (4) 


Concave 


m = 1 


Equilateral Triangle with Interior Vortex (8) 




< m < 1 


Kite 34 (8) 
Asymmetric (8) 




-1/2 < m < 


Kite 19 (4) 


Collinear 


m = 1 


Symmetric (12) 




< m < 1 


Symmetric (4) 
Asymmetric (8) 




-1 < m < 


Symmetric (2) 




-1/2 < m < 


Asymmetric (4) 



Table 1: The number of relative equilibria equivalence classes for the four- vortex problem 
with vortex strengths Ti = T 2 = 1 and T 3 = T 4 = m, in terms of m and the type 
of configuration. The special value m* —0.5951 is the only real root of the cubic 
9m 3 + 3m 2 + 7m + 5. Kite^ refers to a kite configuration with vortices i and j on the axis 
of symmetry. 

masses are assumed to be opposite each other in a convex central configuration [21 [23], is 
still open for the case when equal masses are assumed to be adjacent. In contrast, in this 
paper, we are able to prove that symmetry is required in the case of four vortices with 
two equal pairs of vorticities. In particular, we show that any convex relative equilibrium 
with m > 0, and any concave solution with m < 0, must have a line of symmetry. For 
the convex case, the symmetric solutions are a rhombus and an isosceles trapezoid. In the 
concave case, the symmetric solution is an isosceles triangle with an interior vortex on the 
axis of symmetry. 

A precise count on the number and type of solutions as a function of the parameter 
m is given in Table [TJ A configuration is called a kite if two vortices are on an axis of 
symmetry and the other two vortices are symmetrically located with respect to this axis. 
Kite configurations may either be concave or convex. When counting solutions, we use 
the standard convention from celestial mechanics that solutions which are identical under 
scaling or rotation are considered equivalent. Note that two solutions identical under a 
reflection are counted separately. The full set of solutions for m = 2/5, m = —1/5 and 
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m = —7/10 (excluding any strictly planar configurations identical under a reflection) are 
shown in Figure [TJ 



ra=2/5 







m =— 1/5 




ra = -7/10 




Figure 1: The full set of solutions for m = 2/5, m 
Y\ = T 2 = 1 are denoted by red disks and vortices T 3 



— 1/5 and m = —7/10. Vortices 
r4 = m by green ones. 



When m = 1 (all vortex strengths equal), there are 26 solutions, all symmetric. Al- 
lowing for relabeling, there are only three geometrically distinct configurations: a square, 
an equilateral triangle with a vortex at the center, and a collinear configuration. This 
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is different than the Newtonian case, where an additional symmetric, concave solution 
exists, consisting of an isosceles triangle with an interior body on the axis of symmetry [T] . 
Part of the reason for the contrast is that the equilateral triangle with central vortex is 
degenerate when all vortices have the same strength [18j [22J . An interesting bifurcation 
occurs as m decreases through m = 1, where the equilateral triangle solution splits into 
four different solutions. If vortex 3 or 4 is at the center of the triangle, then the solu- 
tion for m = 1 bifurcates into two different isosceles triangles with the interior vortex 
on the line of symmetry. In this case, the larger-strength vortices are on the base of the 
triangle. If vortex 1 or 2 is at the center of the triangle, the solution branches into two 
asymmetric concave configurations that are identical under a reflection. In this case, the 
weaker-strength vortices are on the base of the triangle. Thus, the number of solutions 
increases from 26 to 34 for the case < m < 1. This concurs with Palmore's result in [22] . 
which specifies a lower bound of 34 non-degenerate solutions for the four-vortex problem 
with positive vorticities. 

As m flips sign, there are two interesting bifurcation values at m = —1/2 and m = 
m* rs —0.5951, the only real root of the cubic 9m 3 + 3m 2 + 7m + 5. There are 26 solutions 
for —1/2 < m < 0, including a family of convex kite configurations having the negative 
strength vortices on an axis of symmetry and the positive strength vortices symmetrically 
located with respect to this axis. Based on the primary result in [3], which is applicable 
to our problem but only when m > 0, one might expect this convex kite to have two axes 
of symmetry, forming a rhombus. However, as demonstrated in Section 7.2 this is not the 
case. There is also a family of symmetric concave kite configurations having one of the 
positive strength vortices in the interior of an isosceles triangle formed by the outer three. 
As m approaches —1/2, the triangle containing the interior vortex and the base of the outer 
isosceles triangle limits on an equilateral triangle while the fourth vortex (at the apex of 
the outer isosceles triangle) heads off to infinity. As m decreases through —1/2, this family 
bifurcates into a convex configuration, while the previous convex kite family disappears. 
This new family of convex kites, with vortices 1 and 2 on the axis of symmetry, exists only 
for m* < m < — 1/2. In sum, there are 18 solutions for m* < m < —1/2, and 14 solutions 
when m = —1/2 or when — 1 < m < m*. The singular bifurcation at m = —1/2 occurs 
in part because the leading coefficient of a defining polynomial vanishes. This is likely 
a consequence of the fact that the sum of three vorticities vanishes when m = —1/2, a 
particularly troubling case when attempting to prove finiteness for the number of relative 
equilibria in the full four- vortex problem (see [IS])- The bifurcation at m = m* is a 
pitchfork bifurcation, as discussed in Section |7.4 



In the next section we define a relative equilibrium and explain how to use mutual 
distances as variables in the four- vortex problem. In Section [3j we describe the relevant 
algebraic techniques used to analyze and quantify the number of solutions. Section [4] ex- 
amines the interplay between symmetry and equality of vorticities in two special cases: 
equilibria and vanishing total vorticity. Sections [5j [6] and [7] cover the collinear, asymmetric 
and symmetric cases, respectively, in considerable detail. Throughout our research, sym- 
bolic computations (e.g., calculation of Grobner bases) were performed using Sage [27], 
Singular [T2] and Maplel5™. 
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2 Relative Equilibria 



We begin with the equations of motion and the definition of a relative equilibrium for 

the n- vortex problem. Let J = ^ _ , and let V,- denote the two-dimensional partial 

— 1 I) 

gradient with respect to Xj. A system of n planar point vortices with vortex strengths 
Tj 7^ and positions X{ £ M 2 evolves according to 

- / rr 

TiXi = JViH = -J^2 ^^(xi-Xj), \<i<n (1) 

where H = — Yli<j ^i^j log 7 "*?? r ij = \\ x i ~ x j\\i an d the prime on the summation indicates 
omission of the term with j = i. 

A relative equilibrium motion is a solution of Q of the form Xi(t) = c+e~ JA '(xj(0) — c), 
that is, a periodic solution given by a uniform rotation with angular velocity A 7^ around 
some point c £ IR 2 . Such a solution is possible if and only if the initial positions Xi(0) 
satisfy the equations 

1 n r 

for each z £ {1, . . . ,n}. Denote T = Yli^i as the total circulation and assume for the 
moment that r / 0. Multiplying the z-th equation in ^ by Tj and summing over i shows 
that the center of rotation c is equivalent to the center of vorticity, c = h J2i ^i x i- If 
T = 0, then we obtain instead that the moment of vorticity £\ T^Xi must vanish. 

Definition 2.1. A set of initial positions Xj(0) satisfying equation ^ for each i £ 
{l,...,n} is called a central configuration. The corresponding rigid rotation with an- 
gular velocity A 7^ is called a relative equilibrium. We will often use these two terms 
interchangeably. 

Define the moment of inertia I with respect to the center of vorticity as 



2 ^ ] Till 3 '* c|| 



i=l 



I measures the size of the system. We can then rewrite equation ^ as 

VH + XV (I - I Q ) = (3) 

where V = (V\, . . . , V n ) and / = I . Therefore, A can be viewed as a Lagrange multiplier 
and any solution of ^ can be interpreted as a critical point of the Hamiltonian H(x) 
under the condition that / remains constant. Using the homogeneity of the functions H 
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and J, equation ^ implies that the angular velocity A in a relative equilibrium is given 
by 

L 

X = - where L = X> r i ( 4 ) 

i<j 

is the total vortex angular momentum. This implies that in the case where all vortex 
strengths are positive, A > and the relative equilibrium is rotating in the counterclockwise 
direction. It is also important to note that for any relative equilibrium with a particular set 
of vortex strengths Tj, we can scale the vorticities by some common factor v and maintain 
the relative equilibrium, but with a new angular velocity v\. If v < 0, then the sign of A 
flips, as does the direction of rotation. 

Unlike the ra-body problem, it is possible to have equilibria in the n-vortex problem 
(where A = 0). For n = 4, explicit solutions can be derived (see [16]). An analysis of the 



equilibrium solutions in the context of symmetry is presented in Section 4.1 



2.1 Using mutual distances as coordinates 



We now consider the case of n = 4 vortices. Our presentation follows the approach of [21] 
in describing the work of Dziobek [13] for the Newtonian n— body problem. We want to 
express equation ^ in terms of the mutual distance variables r^. Between four vortices 
there are six mutual distances, which are not independent if the vortices are planar. In 
the planar case the distances satisfy the following condition, which can be interpreted as 
the vanishing of the volume of the tetrahedron formed by the four vortices: 






1 


1 


1 


1 


1 











1 


r u 





iy 2 

' 23 


iy 2 

' 24 


1 




'23 





/y> 2 

'34 


1 


r 2 u 


/>i 2 
'24 


,y> 2 

'34 






The matrix in the above determinant is known as the Cayley-Menger matrix. 
Hence, planar central configurations are obtained as critical points of 



f+a(j-/ ; 



(5) 



Using the homogeneity of H, I and ecM, the value of A in this setup is identical to one given 
in equation Q. To find VecM restricted to planar configurations, we use the following 
important formula 

de C M 



where Ai is the oriented area of the triangle Tj whose vertices are all the vortices except 
for the i-th body. Setting the gradient of ^ equal to zero yields the equations 

de C M 



OH x dl 
+ A 



dr 2 - 



32 dr\ 



0. 



7 



If r 7^ 0, then / can be written in terms of the mutual distances as 

i = — Vr r r 2 . 

i<j 

Using this, we obtain the following equations for a four-vortex central configuration: 

T l T 3 {r- 2 + \') = aA l A 3 (6) 

where A' = — A/T, o = — 2/z, I = I and ecu = 0. If Y — 0, then a different approach is 
more useful (see [8j [16] for expositions of the equations). We discuss the role of specific 



symmetries when T = in Section 4.2 



Assuming r ^ 0, we group the equations in (J6J) so that when they are multiplied 
together pairwise, their right-hand sides are identical: 

rxr 2 (rr 2 2 + A') = aA x A 2 , T 3 T 4 {r^ + A') = aA 3 A 4 

T.T^ 2 + A') = aA x A z , T 2 T 4 {r 24 2 + A') = aA 2 A 4 (7) 

rir 4 (rr 4 2 + A') = aA x A A , T 2 T 3 (r^ 2 + A') = crA 2 A 3 . 

This yields the well-known Dziobek equations [13J, but for vortices: 

(rr 2 2 + A')^ 2 + A') = (rr 3 2 + A')(r 2 " 4 2 + A') = (r u 2 + A')^ 2 + A'). (8) 

From the different ratios of two masses that can be found from the equations in (J7| we 
obtain the following equations: 



t x a 2 


_ P23 + A' 


_ P24 + A' 


_ P23 


- P24 


t 2 a 1 


Pis + A' 


Pl4 + A' 


Pl3 


- Pl4 


t x a 3 


_ P23 + A' 


_ P34 + A' 


_ P23 


- P34 


T 3 A 1 


P12 + A' 


Pl4 + A' 


Pl2 


- Pl4 


t x a 4 


_ P24 + A' 


_ P34 + A' 


_ P24 


- P34 


r 4 A x 


Pl2 + A' 


Pl3 + A' 


Pl2 


- Pl3 


T 2 A 3 


_ P13 + A' 


_ P34 + A' 


_ Pl3 


- P34 


T 3 A 2 


Pl2 + A' 


P24 + A' 


Pl2 


- P24 


T 2 A 4 


_ Pl4 + A' 


_ P34 + A' 


_ Pl4 


- P34 


T 4 A 2 


Pl2 + A' 


P23 + A' 


Pl2 


- P23 


t 3 a 4 


_ Pl4 + A' 


_ P24 + A' 


_ Pl4 


- P24 


r 4 A 3 


Pl3 + A' 


P23 + A' 


Pl3 


- P23 



(9) 



where pij = r. 
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Eliminating A' from equation ([8]) and factoring yields the important relation 

(r 2 13 ~ r 2 2 )(r 2 3 - r 2 4 )(r 2 4 - r 2 4 ) = (r 2 2 - r 2 4 )(r 2 4 - r 2 4 )(r 2 3 - r 2 3 ). (10) 

Assuming the six mutual distances determine an actual configuration in the plane, this 
equation is necessary and sufficient for the existence of a four-vortex relative equilibrium. 
The corresponding vortex strengths are then found from the equations in 
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Relationships between the lengths of the sides in a four-vortex central configuration 
follow from the equations in ^ and the signs of the oriented areas A> If the configuration 
of vortices is concave, precisely three of the oriented areas have the same sign. In the 
convex case, two oriented positive and two are negative. 

The sign of A' = — A/T for a given relative equilibrium is of some interest. Note that 
scaling the vorticities by any v G IR — {0} does not change the value of A'. If all the 
vorticities have the same sign, then A' < is assured. However, if the vorticities have 
different signs, it is possible that A' could become positive. 

When A' > 0, the equations in (9 ) imply that ^-r- > for any choice of indices i and j. 
Taking T 1 = T 2 = 1 and T 3 = T 4 = m with m < 0, we see that the only possible solutions 
have Ai,A 2 > and A3, A < 0, or A, A < and A, A > 0. The configuration must 
be convex with vortices 1 and 2 on one diagonal and vortices 3 and 4 on the other. We 



show in Section 7.3 that the configuration must have at least one axis of symmetry. There 



exists a fami ly o f convex kite configurations for m* < m < — 1/2 and a family of rhombi 



see Section 7.4) for — 1 < m < — 2 + \^3. These are the only possible solutions to our 



problem having A' > 0. 

Proposition 2.2. Suppose we have a concave central configuration with four vortices and 
r ^ 0. 

1. If all the vorticities are positive (negative), then all exterior sides are longer than the 
interior ones and \/\J—X' is less than the lengths of all exterior sides and greater than 
the lengths of all interior sides. 

2. If two of the vorticities are positive and two are negative, then the exterior sides con- 
necting vortices with vorticities of opposite sign and the interior side connecting vortices 
with vorticities of the same sign have lengths greater than 1 / y/—X', while the remaining 
sides have length less than l/y/—X'. 

Proof. (1) Let be the interior vorticity and suppose that A < 0. Then we have 
Ai, A2, A 3 > 0. Assume that P34 > —A'. Then, using the equations in Q, we obtain 

Pl2,Pl3,P23 < "A' < Pl4,P24,P34, 



or 

ri2,ri3,r 2 3 > 1/v^A 7 > r u ,r 2 4,r u , 

namely all the exterior edges are longer than the interior ones. On the other hand, if we had 
assumed P34 < —A', then the inequalities above would be reversed and the configuration 
could not be realized geometrically since the interior sides cannot all be longer than the 
exterior sides. 

(2) Let T 4 be the interior vorticity. Let r 3 ,r 4 < 0, and ri,r 2 > 0. Furthermore, 
suppose that At < 0. Then we have Ai,A 2 ,A 3 > 0. If p 34 < —A', then, using the 
equations in Q, we obtain 

Pl3,P23,P34 < -A' < Pl2,Pl4,P24, 
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or 

^13,^23,^34 > l/v /r A 7 > ri2,r 14 ,r 2 4. 

If we had assumed P34 > — A', then the inequalities above would be reversed and the con- 
figuration could not be realized geometrically. We show this with a proof by contradiction. 
Assume pu > — A'. Then \A±\ < \As\ since the two triangles have one side in common and 
^24,^14 > ^23,^13- But this is absurd since T4 contains T3. 

Let Ti be the interior vorticity. Let T3,r4 < 0, and Ti,r2 > 0. Furthermore, suppose 
that Ai < 0. Then we have A 2 , A 3 , A 4 > 0. If p 34 > —A' then, using the equations in (|9]), 
we obtain 

Pl2,P23,P24 < -A' < Pl3,Pl4,P34, 

or 

^12,^3,^24 > 1/v^A 7 > ri 3 ,ri4,r 3 4. 

If we had assumed p 34 < —A' then the inequalities above would be reversed, and the 
configurations could not be realized geometrically (the proof is similar to the argument in 
the previous paragraph). □ 

Before we state an analogous theorem for the convex case we recall the following useful 
geometric lemma: 

Lemma 2.3. The combined length of the diagonals of a convex quadrilateral is greater 
than the combined length of any pair of its opposite sides. 

Proof. Order the vortices counterclockwise. Let o be the position vector describing the 
intersection of the diagonals, and let r io = \\xi — o\\. Applying the triangle inequality to 
the triangle of vertices x±, x% and o and to the one of vertices 23, £4 and o yields 

r%o + r 2o > r i2 and r 3o + r 4o > r 34 . 

Adding these two inequalities together, we obtain r 13 + r 2 4 > 7" 12 + r3 4 . A similar reasoning 
can be applied to the remaining two triangles. □ 

Proposition 2.4. Suppose we have a convex central configuration with four vortices, with 
r^O and A' < 0. 

1. If all the vorticities are positive (negative), then all exterior sides are shorter than the 
diagonals. Furthermore, the lengths of all exterior sides are less than 1/v — A' and the 
lengths of all the diagonals are greater than 1 / y/—X'. The shortest and longest exterior 
sides have to face each other. 

2. If two of the vorticities are positive and adjacent, and two remaining ones are negative, 
then the exterior sides connecting vortices with vorticities of opposite sign have length 
less than 1/v — A 7 . All the other sides have length greater than l/y 7 — A 7 - 

3. If two of the vorticities are positive and opposite, and two are negative, then either 
all the sides have length less than 1/V— A' , or all the sides have length greater than 
l/v^. 
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Proof. (1) The proof is analogous to the one for the Newtonian four-body problem (see 

(2) Order the vortices counterclockwise and let r 1; T 2 > and T 3 , T 4 < 0. The sign of 
the four areas are then A\ > 0, A 3 > and A 2 < 0, < 0. If we assume that p 34 > —A', 
then we find from the equations in ^ that 

Pl2,Pl3,P24,P34 < -A' < Pl4,P23, 

or 



r 12, r '13, r 24 ,r 34 > 1/V-A' > ri 4 ,r 23 . 
On the other hand, if we had assumed that p 34 < —A', then the inequalities above would 



be reversed and, by Lemma 2.3, the configuration could not be realized geometrically. 

(3) Let r 1; r 2 > and T 3 , T 4 < 0, and assume the vortices of the same sign are opposite 
one another. The sign of the four areas are then A\ > 0, A 2 > and A3 < 0, A 4 < 0. If 



we assume that p 34 > — A', then we find from the equations in (|9j) that 

P12, P13, P24, P34, Pu, P23 > —A', 

or 

^12,^13,^24, ^34,^14,^23 < V V^X ■ 

If p 34 < —A', then the inequalities are reversed. □ 
2.2 Symmetric configurations 



One immediate consequence of the Dziobek equations (|8j) and of equation (10) is that if 
two mutual distances containing a common vortex are equal (e.g., r 12 = r 13 ), then the 
same equality of distances is true for the excluded vortex (here, r 24 = r 34 ). Specifically, if 
i, j, k, I are distinct indices, then we have 

Tij = r ik if and only if r tj = r lk . (11) 

This relation is independent of the vortex strengths although T,- = T k necessarily must 



follow. Any configuration satisfying equation (11) has an axis of symmetry containing 
vortices i and /, forming a kite configuration. It may be either convex or concave, but it 
cannot contain three vortices on a common line due to the equations in ([7]). Unlike the 
Newtonian four-body problem, since there is no restriction here on the signs of the vortex 
strengths, it follows that any kite configuration has a corresponding set of vorticities that 
make it a relative equilibrium. Assuming the configuration is not an equilateral triangle 
with a vortex at the center, the vorticities are unique up to a common scaling factor and 
can be determined by the equations in 

Next we consider the case where two mutual distances without a common index are 
equal. In other words, suppose one of the following three equations holds: 

ri2 = r 3A , r 13 = r 24 , r u = r 23 - (12) 
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Then, it does not necessarily follow that another pair of mutual distances must be equal. 
However, if different pairs of vortices are assumed to be of equal strength, then we can 



conclude an additional symmetry. This fact will be important in Section 6T when verifying 
that symmetry is required in certain cases. 

Lemma 2.5. (Symmetry Lemma) Suppose that we have a strictly planar four-vortex 
relative equilibrium with Ti = T 2 , T 3 = T 4 and r ^ 0. Then, 

ria = r M if and only if r 14 = r 23 . (13) 



If either equation in (15) holds, the configuration is convex and has either one or two axes 
of symmetry. In this case, the configuration is either an isosceles trapezoid with vortices 
1 and 2 on one base, and 5 and 4 on the other, or it is a rhombus with vortices 1 and 2 
opposite each other. 

Proof. Suppose that r 13 = r 24 . Then, by one of the equations in Q, r 1 r 3 (r^ 3 2 + A') = 
r 2 r 4 (r^ 4 2 + A') implies that 

A X A 3 = A 2 A A . (14) 

We also have that 

A 1 + A 2 + A 3 + A 4 = (15) 



since the A^s are oriented areas. Solving equation (15) for A 4 and substituting into 



equation (14) yields the relation 

{A l + A 2 ){A 2 + A 3 ) = 0. 



There are two possibilities. First, suppose that A% = —A 2 . Then equation (15) imme- 
diately implies A 4 = — A 3 . The configuration is convex due to the signs of the A^s and 
must have the side containing vortices 1 and 2 parallel to the side containing vortices 3 
and 4. Then, by another equation in ([7]), aA 1 A 4: = aA 2 A 3 implies that 

rir 4 (rr 4 2 + A') = r 2 r 3 (r 2 - 3 2 + A'). 

Since I\ = T 2 and T 3 = T 4 , it follows that ri 4 = r 23 . The configuration must be an isosceles 
trapezoid with congruent legs and diagonals, and base lengths given by r\ 2 and r 34 . 

Next, suppose that A 3 = —A 2 . Then, v4 4 = —A\ immediately follows from equa- 



tion (15). The configuration is convex and must have the side containing vortices 1 and 4 
parallel to the side containing vortices 2 and 3. Since r 13 = r 24 , the configuration is either 
an isosceles trapezoid with r\ 2 = r 34 or is a parallelogram with ru = r 23 . For the isosceles 
trapezoid, the lengths r i3 and r 24 can correspond to either a pair of congruent diagonals 
or a pair of congruent legs, depending on the ordering of the vortices. However, for the 
parallelogram, these lengths must correspond to a pair of opposite sides as the diagonals 
do not have to be congruent. 

It turns out that, in this case, the isosceles trapezoid is actually a square. To see this, 
by an equation in ([7]), o~A\A 2 = oA 3 A± implies that 

r 2 (rr 2 2 + A') = r 2 (r 3 - 4 2 + A'). 
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Thus, if r 12 = r 34 , then Ti = T 3 and all vortex strengths are equal (the case I\ = — T 3 
is excluded since f ^ 0). In this case, it is straight-forward to show that the isosceles 



trapezoid reduces to a square (see Section 7.1) and thus r i4 = r 23 . 



In the case of the parallelogram, we have A\ = A 2 = —A3 = — A4. By the equations 
in (J7|), this implies that ri 3 = ri 4 = r 23 = r 24 and the configuration is a rhombus with 
vortices 1 and 2 opposite each other. This proves the forward implication. 

The proof in the reverse direction is similar. If r 14 = r 23 , then we derive A\A± = A 2 A 3 
from an equation in ([7]). Taken with equation (15), this yields 



(A 1 + A 2 )(A 1 + A 3 ) = 0. 

As before, the case A\ = —A 2 leads to an isosceles trapezoid with r 13 = r 24 and base 
lengths given by r\ 2 and r 34 . The case A\ = — v4 3 leads to either the square or a rhombus 
configuration with vortices 1 and 2 across from each other. In either configuration we 
deduce that r 13 = r 24 . This completes the proof. □ 



Remark. 1. Similar results exist if different pairs of vortices are assumed to be equal. 
For example, if I\ = T 3 and T 2 = T 4 , then r 12 = r 34 if and only if r 14 = r 23 . 

2. The result is also valid in the Newtonian four-body problem (and for other potentials 
of the same form) since it only depends on the geometry of the configuration and the 
inherent structure of the equations in ([7]). 

2.3 The Albouy-Chenciner equations 

For the remainder of the paper (excluding the special cases discussed in Section 111), we 



will assume the equality of vortex strengths specified in Lemma 2.5 Specifically, we set 
Ti = T 2 = 1 and T 3 = T 4 = m, treating m as a real parameter. Without loss of generality, 
we restrict to the case where m G (—1, 1]. The choice m = —1 implies that T = 0, a 
special case examined in Section |4.2[ 

When r 0, the equations for a relative equilibrium can be written in polynomial 
form as 

n 

fa = T kl S ik(r 2 jk - rj k - rj) + S jk {r 2 ik - r) k - r\)\ = 0, 

k=l 

where 1 < i < j < n and the Sy are given by 

Sij = -j + X' (i^ j), s u = 0. 

These very useful equations are due to Albouy and Chenciner [2] (see also [15] for a 
nice derivation). They form a polynomial system in the ry variables after clearing the 
denominators in the terms. 
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Since any relative equilibrium may be rescaled, we will impose the normalization 
A' = — 1 unless otherwise stated. This usually can be assumed without loss of gener- 
ality. However, as explained in Section 2.1\ for m < the normalization A' = 1 also needs 



to be considered. This case is discussed in Sections |7.3| and |7.4| We denote the complete 
set of polynomial equations determined by = as T . 

From the Albouy-Chenciner equations we can derive a more restrictive set of equations, 
namely 

n 

9a = Yl T k S ik(r 2 jk - r] k - rj) = 0. 

k=l 

Since ^ gji, these give 12 distinct equations. We denote the complete set of polynomial 
equations determined by g^ = as Q, and we will refer to them as the unsymmetrized 
Albouy-Chenciner equations. 

The solutions of the Albouy-Chenciner equation give configurations of all dimensions, 
but, in the four- vortex problem, they can be specialized to the strictly planar case by 
adding the three Dziobek equations in dsb - Introducing the variables = r?-, these 
equations can be written as 

h m = {s$ + \')(s- k ? + A') - fo 1 + A')^ 1 + A') = 0, 

where i,j, k and I are all distinct indeces. We denote the set of Dziobek equations (with 
denominators cleared) as %. 



3 Algebraic Techniques 

In this section we briefly describe three of our main algebraic techniques for analyzing 
solutions to our problem: elimination theory using Grobner bases, a useful lemma to 
distinguish when the roots of a quartic are real or complex, and Mobius transformations. 

3.1 Grobner bases and elimination theory 

We mention briefly some elements from elimination theory and the theory of Grobner 
bases that will prove useful in our analysis. For a more detailed exposition see [TT] . 

Let K be a field and consider the polynomial ring K[ ] of polynomials in n 

variables over K. Let fi, ■ ■ ■ fi be / polynomials in K[x\, . . . ,x n ] and consider the ideal 
I = (fx, . . . , fi) generated by these polynomials. Denote V(J) as the affine variety of /. 

Definition 3.1. An admissible order > on K[x±, . . . ,x n ] is called a k— elimination order 
if 

x l • • • ^ ' ' ' n 

when a io > for some zo G {1, . . . , k}. 

A lexicographical (lex) order is an example of k— elimination order for all k. 
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Definition 3.2. The k-th elimination ideal Ik is the ideal of K[xk+i, ■ ■ ■ ,x n ] defined by 

h = ir\ K[x k+1 , ...,x n ] 

Grobner bases provide a systematic way of finding elements of I k using the proper term 
ordering. 

Theorem 3.3. (The Elimination Theorem) Let I be an ideal of K[ Ob \ . • • • . Ob fl 1 and let 

G be a Grobner basis of I with respect to a k— elimination order for k where < k < n. 
Then the set 

G k = Gn K[x k +i, ■ ■ ■ ,x n ] 
is a Grobner basis of the k-th elimination ideal Ik- 

Grobner bases also provide a method of determining when an element of V(/fc) (a 
partial solution) can be extended to a full solution in V(J). This can be achieved by 
repeatedly applying the following theorem. 

Theorem 3.4. (The Extension Theorem) Let K be an algebraically closed field and 
let I be some ideal in K[xi, . . . ,x n ]. Let Gk-i be a lex Grobner basis for the elimination 
ideal Ik-i and write each polynomial in Gk-i as 

gi = hi(xk+i, ■ ■ ■ Xrij^k* + t erms where x k has degree < JVj, 

where Ni > and hi is nonzero. Suppose that (a,k+i, • • • , a n ) is a partial solution in V(/fc) 
and that hi(ak+i, ■ ■ ■ , a n ) ^ for some index i. Then there exists G K such that 
(a k ,a k+ i, ...,a n )e V(4_i). 

3.2 A useful lemma for quart ic polynomials 

The analysis of the collinear well as some strictly planar cases, frequently involves 

solving a quartic equation whose coefficients are polynomials in m. We state here some 
useful results about quartic equations. Consider the general quartic polynomial ((x) = 
ax A + bx 3 + cx 2 + dx + e, with coefficients in R and a / 0. We first remove the cubic term 
of ( by the change of variables x = y — ^. This produces the polynomial arj(y), where 77 
is the shifted quartic rj(y) = y A +py 2 + qy + r. The discriminant A of any polynomial is a 
positive constant times the square of the product of all possible differences of roots. The 
discriminant of £ is equivalent to a 6 times the discriminant of r]. For a general quartic, it 
is straight-forward to check that if A > 0, then the roots are either all real or all complex 
(two pairs of complex conjugates). If A < 0, then there are two real roots and two complex 
roots. The roots are repeated if and only if A = 0. 

Let yi,y2, 2/3 and 7/4 be the four roots of rj(y). By construction, 2/1+2/2 + 2/3 + 2/4 = 0. 
It follows that 

Zi = -(2/1 + 2/2) 2 , z 2 = -(2/1 + 2/3) 2 , z 3 = -(y 1 + 2/ 4 ) 2 
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are the roots of the resolvent cubic ^(z) = z 3 — 2pz 2 + (p 2 —Ar)z+q 2 . The discriminant of £ is 
equivalent to the discriminant of 77. When the discriminant is positive, the resolvent cubic 
is particularly useful for determining whether the roots of rj(y) are all real or all complex. 
Specifically, if the four roots of T](y) are real, then zi,z 2 , and 23 must be negative and real 
(or if q = 0, then Zi = for precisely one i while the remaining Zj's are negative). On the 
other hand, if the four roots of r](jj) are complex, then one of the roots of £(z) is less than 
or equal to zero, but two of the roots of £(z) must be positive and real. These facts can 
easily be translated into conditions on the coefficients of the resolvent cubic. 

Lemma 3.5. Suppose that the discriminant of a quartic polynomial ((x) is positive and 
let r)(y) = y A +py 2 + qy + r be the shifted quartic related to (. Then, the four distinct roots 
of ( are real if and only if p < and p 2 — 4r > 0. 

Proof. Since the discriminant is positive, the roots of ( are distinct and either all real or 
all complex. First suppose that the four roots of ((x) are real. Since rj(y) is obtained 
through a simple translation, it follows that the four roots of rj(y) are distinct and real. 
Then, the three real roots z\, 2:2,23 of the resolvent cubic £(z) are negative with possibly 
one root equal to zero (if q = 0). This in turn implies that 

zi + z 2 + Z3 = 2p < and 
Z\Z 2 + ZxZz + z 2 z 3 = p 2 - 4r > 0. 

In the other direction, suppose that p < and p 2 — 4r > 0. This implies that all the 
coefficients of £(z) are positive, and by Descartes' rule of signs, £(z) has no positive, real 
roots. It follows that the four roots of rj(y) cannot be complex. This shows that the roots 
of ((x) are real. □ 



3.3 Mobius transformations 

Mobius transformations are a key tool for isolating roots of polynomials since they allow 
for a dramatic reduction on the number of variations of signs in the coefficients of a 
polynomial [5]. Let V(x±, . . . , x n ) be a multivariate polynomial in n variables. We use 
Mobius transformations of the form 

where kf^ < and i — 1, . . . n, to obtain changes of variables for V(x\, . . . , x n ). The 
numerator of the rational function thus obtained is a multivariate polynomial in the vari- 
ables yi, . . . , y n that restricted to the set [0, 00) x . . . x [0, 00) has the same number of 
roots as V(xi, . . . ,x n ) restricted to [k^, k^) x . . . x [kn \ kn^). 

Mobius transformations can be used to determine partitions of the space into blocks 
so that, in each block, the coefficients of the polynomials under consideration have a 
very simple behavior with respect to variations of signs. Therefore, we can use such 
transformations, together with Descartes' rule of signs, to determine whether a polynomial 
has a root or does not have a root in a given region of space. 
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4 Special Cases 



In this section, we use some results from [16J to examine two special cases: equilibrium 
solutions and vanishing total vorticity (T = 0). 



4.1 Symmetric equilibria 

Equilibria are solutions to equation ^ with A = 0. If we choose coordinates so that the 
fourth vortex is at the origin, and the third vortex is at (1, 0), then the other two vortices 
in a four-vortex equilibrium must be located at: 

* = 2(r, + r, + ro fo + r„±^r,) ^ 
^ 2 ( r 1 + r 3 + r 4 ) fo + r^r.) 

where the sign chosen for x 2 is the opposite of that for x\ (see [E] for details). 

According to equation Q, for an equilibrium to exist the vorticities must satisfy L = 0. 
This is impossible if all the vorticities are equal. If three vorticities are equal, for example 
with Ti = T 2 = T 3 = m, then T 4 = —m. We can rescale so that m = 1 without loss of 



generality. Using (16) it is easy to see that the only equilibria are equilateral triangles 
with the opposing vortex located in the center. 

Now consider equilibria when Ti = T 2 = 1 and r 3 = r 4 = m. Then, L = implies 
m = — 2 ± v^3- Since we are restricting to m E (—1, 1], we have m = — 2 + ~ —0.2679. 
In this case, the equilibria are rhombi with vortices 1 and 2 opposite each other, and 
— = 2 — V3 = —m. These rhombi are members of one of the families of rhombi described 
in Section \Tm 



Suppose we instead specify the symmetry of the configuration. From (16) it is im- 
mediate that there are no collinear equilibria for nonzero vorticities. There are also no 
isosceles trapezoid equilibria. To see this, consider a trapezoid with r 14 = r 2 3 and r 13 = r 2 4. 



Then (16) implies that 

r 2 -ri 



r 2 + r 3 + r 4 r 1 + r 3 + r 4 ' 

but there are no real nonzero vorticities satisfying this equation and L = 0. 

Finally, there is the case in which the configuration is a kite (either concave or convex) . 
We choose vortices 1 and 2 to be on an axis of symmetry, so ri 3 = ru and r 23 = r 24 . Then 



equation (16) implies that T 3 = T 4 . If we choose any T 2 ^ — 2r 4 with 



2rxr 4 + r 2 



4 



v l + 2r 4 



then L = is satisfied and there is a kite equilibrium given by (16). If we fix T 4 = 1, 
then as T 2 — > —2, Ti — )• oo and the configuration of vortices 2, 3, and 4 approaches an 
equilateral triangle while the ordinate of vortex 1 heads off to ±oo. 
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4.2 Vanishing total vorticity 



When the total vorticity T = 0, the analysis of stationary vortex configurations (i.e., those 
that do not change their shape) requires equations adapted to this special case. For four 
vorticities we simply apply some of the results from [16J. 

If three of the vorticities are equal and T = 0, we can consider without loss of generality 
the case Fx = T 2 = T 3 = 1 and T 4 = —3. It is interesting that there are six asymmetric 
configurations which rigidly translate (three pairs of configurations, within a pair the con- 
figurations are reflections of one another), but no symmetric rigidly translating solutions. 
These configurations can be obtained directly from the equations in [16J. 

For relative equilibria we use the equations 

Si — 02 — S3 — 64 = So 

and 

111111 

— + — = — + — = — + —, (17) 

Sl2 S34 S13 S 2 4 S14 S 2 3 

where 

Si = TjSij + T k s ik + Tis u , {i,j,k,l} = {1,2,3,4} 
and So is an auxiliary variable. Recall that = rh. We clear denominators in the 



equations from (17) to get a polynomial system. Using a Grobner basis to eliminate sq, 
we find that for Fx = T 2 = T 3 = 1 and T 4 = —3 there are two types of symmetric 
relative equilibria. The first is the equilateral triangle with vortex 4 at its center. The 
second type is a concave kite with the three equal vorticities on the exterior isosceles 
triangle. If we scale the exterior triangle so that its longest side is length 1, then the base 



is length y y3 — 5, and the other sides of the interior triangle containing vortex 4 are 
length y/-l + y/3. 

If two pairs of vorticities are equal and T = 0, then we have Ti = T 2 = 1 and 
T3 = T 4 = — 1. In this case there are no rigidly translating solutions. This is somewhat 
surprising since each pair of opposing vortices would rigidly translate if unperturbed. Using 
the same equations as we did for the case where three vorticities are equal, we find that 



there are two relative equilibria, each o: 
families of rhombi described in Section 



which forms a rhombus. These appear in the two 

with — = 3 ± 2V2. 
S12 



7.4 



5 Collinear Relative Equilibria 

Collinear relative equilibria of the four- vortex problem can be studied directly from equa- 
tion g since in this case it reduces to 

n p 

-\( Xi -c) = Y^ — for each i G {1,2,3,4}, (18) 

where c,Xi G K rather than IR 2 . Clearing denominators from these equations yields a 
polynomial system. Rather than fix A or c, we use the homogeneity and translation 
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invariance of the system and set x 3 = — 1 and x 4 = 1. As specified earlier, we set 
r\ = r 2 = 1 and r 3 = T 4 = m G (—1, 1], treating m as a parameter. 



5.1 Symmetric solutions 

Given our setup, symmetric configurations correspond to solutions where either x% = —x 2 
or 7*12 — 2. Both cases are simple enough to analyze using Grobner bases. 

Case 1: x\ = —X2 

In this case the center of vorticity c is at the origin. The solutions are the roots of an even 
quartic polynomial in x%, given by mxf — 5(m + l)x\ + 1. The discriminant of this quartic 
changes sign at m = and there are four real roots for m > but only two when m < 0. 
The result is that for any m G (—1,1], there is a collinear relative equilibrium with 



/ 5 m + 5 - V25m 2 + 46m + 25 

Xl = - X2 = ±y (19) 

(with Xi = —X2 = ±l/Vo when m = 0). This solution has vortices X\ and x 2 symmetri- 
cally located between vortices £3 and X4, and as m — >■ — 1 + , the inner vortices approach 
the outer vortices (with collision at m = —1.) For m G (0, 1], \x\\ decreases monotonically 
in m from I/a/5 ~ 0.4472 to \/3 — \/2 w 0.3178, while for m G (—1,0), |a;i| decreases 
monotonically in m from 1 to 1/ yfb. 

There is an additional collinear relative equilibrium if m G (0, 1] given by 



X\ = —x 2 = ±< 



'5m + 5 + v / 25m 2 + 46m + 25 
2m 



In this case, the vortices X\ and x 2 are symmetrically located outside vortices x% and x±, 
and their positions approach ±00 as m — » + . The value of \x\\ decreases monotonically 
m m from 00 to V3 + ^2 w 3.1462. 

Case 2: r i2 = 2 

In this case the center of vorticity c is not necessarily located at the origin; however, the 
configuration will be symmetric about some fixed point equidistant from both the inner 
and outer pairs of vortices. Since \x\ — x 2 \ =2 and x\ = —x 2 together imply \x 2 \ = 1 
(collision), the two cases are distinct. If m ^ 1, a solution for case 2 will have the inner and 
outer pair of vortices having different circulations. It turns out that this case is impossible 
when m/1. 



To see this, we consider the polynomials obtained from system (18) along with x±— x 2 — 2 
and u{x\ + x 2 ) — 1 (to eliminate solutions from Case 1). Computing a lex Grobner basis 
for this set of polynomials quickly yields m = 1. The same result is obtained when using 
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Xi — X2 + 2. Moreover, in each computation, a fourth-degree polynomial in x 2 is obtained 
that provides the exact solution for the special case m—1. 

In sum, there are no solutions with r X2 = 2 unless m = 1. When m = 1, there are 8 
solutions with r X2 = 2 given by (xi,x 2 ) = 

(1±V2±V6,-1±V2±V6) and (-1 ± y/2 ± y/E, 1 ± V2 ± V6), (20) 

where the signs in a particular ordered pair are chosen to have the same pattern in each 
coordinate (e.g., (1 + — \/6, — 1 + \[2 — \/6)). It is straight-forward to check that each 
of these 8 solutions is a scaled version of a solution found in Case 1. Specifically, for each 



of the solutions listed in (20) as well as the solutions obtained in Case 1 for m—1, there 



is a rescaling, translation and relabeling of the vortices that maps the solution onto 

x 3 = -l, x 2 = -V3 + V2, x x = V3-V2, 2:4 = 1, (21) 



a solution obtained from equation (19) when m—1. We note that solution (21) can be 



rescaled to coincide with the roots of the Hermite polynomial H4, as expected (see [1]). 
5.2 Asymmetric solutions 

To locate any asymmetric solutions, we introduce the variables u and v along with the 
equations 

u{x\ + X2) — 1 and v{x\ — x 2 ) — 1. 



Adding these two equations to the original polynomial system obtained from (18), we 
compute a Grobner basis G co \ with respect to the lex order where c>A>w>i;>£i> 
x 2 > m - This basis has 15 elements, the first of which is an even, 8th-degree polynomial in 
X2 with coefficients in m, and a basis for the elimination ideal I5 = 7flC[x2, m\. Introducing 
the variable w = x\, this polynomial is given by 

((w) = m 2 (m + 2)(1 + 2m) 2 w 4 — 4m(15m 4 + 61m 3 + 91m 2 + 61m + 15) ^ 3 
+ (300m 5 + 1508m 4 + 2910m 3 + 2696m 2 + 1188m + 200) w 2 
-4(5m + 4) (25m 4 + 127m 3 + 231m 2 + 175m + 45) w + (m + 2) 3 . 

If we eliminate x 2 instead of x%, the same polynomial is obtained with w = x\. The 
discriminant of ( is 

1048576(m + 2) 2 (m + 1) 6 (1 + 2m)(25m 2 + 58m + 25) 3 m 2 (g n (m)) 2 

where q u {m) = 2m 5 — 16m 4 — 96m 3 — 162m 2 — 108m — 25. For m e (0, 1], the discriminant is 
strictly positive. However, on the interval [—1, 0], the discriminant vanishes at six different 
m-values. In increasing order, these values are 

-1, m w -0.6833, m x w -0.6066, m 2 w -0.5721, -1/2, 0. 

The values mo and m\ are roots of the quintic polynomial q u , and m 2 = (—29 + 6v6)/25 
is the largest root of the quadratic 25m 2 + 58m + 25. For m = — 1, ( has a repeated root 
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at 1 of multiplicity four. At m — —1/2, ( becomes a cubic polynomial with roots —3 and 
1 (multiplicity 2). Thus, the values m = — 1 and m = —1/2 have no physical solutions 
(just collisions). For m = 0, the quartic reduces to a quadratic, and we have two possible 
solutions at w = (9 ± Ay/5)/ 5 which correspond to limiting configurations as m — > + . 
The other three values where the discriminant vanishes also have no physical solutions, as 
we show below. 

Lemma 5.1. The quartic polynomial C( w ) has precisely four positive real roots for m G 
(0,1] and two positive real roots form G (—1/2,0). Form G (—1,-1/2), there are no 
positive roots except when m = m , where there is one positive, repeated root of multiplicity 
two. 



Proof. Using Lemma 3.5, let rj be the shifted quartic of ((w). From rj, the key quantities 
p and p 2 — 4r are found to be 

^ _ 2(25m 2 + 58m + 25) (3m 4 + 14m 3 + 45m 2 + 54m + 19) (m + l) 2 



and 



p 2 — 4r 



m 2 (m + 2) 2 (1 + 2m) 4 
16(25m 2 + 58m + 25) (m + l) 2 t(m) 



m 4 (m + 2) 4 (l + 2m) 8 
where 

t(m) = 75m 12 + 896m 11 + 5528m 10 + 23492m 9 + 77272m 8 + 197816m 7 + 376194m 6 
+509968m 5 + 478976m 4 + 302388m 3 + 121964m 2 + 28320m + 2875. 



For m G (0, 1], since the discriminant is positive, p < and p — 4r > 0, Lemma 3.5 
applies and the roots of C( w ) are all real. The fact that they are all positive follows 
from Descartes' rule of signs since ({w) has four sign changes while ((— w) has none. For 
m G (—1/2,0), the discriminant remains positive and using Sturm's theorem, one can 
show that p < and p 2 — 4r > both continue to hold. It follows that ((w) has four 
real roots in this case as well. However, two of the roots are positive and two are negative 
since both ((w) and £(— w) each have two sign changes. 

For the case m G (m2, —1/2), the discriminant is negative and thus there are precisely 
two real roots. To show that the real roots are both negative, we note that at m = mi, the 
quartic ( has a repeated root of multiplicity four at w 2 = — (11 + 4 a/6 )/5 < 0. We then 
can check that Ci w 2) as a function of m is strictly negative for m G (m 2 , —1/2). Since the 
leading coefficient of £ is positive and since ((0) = (m + 2) 3 > 0, it follows that both real 
roots must be negative in this case. There are no positive roots at m = m 2 since W2 < 
is the only root. 

For the interval m G (— l,m 2 ), the discriminant is positive except when it vanishes at 
mo and m\. However, either p > or p 2 — 4r < (or both) for m- values on this interval. 



Consequently, by Lemma 3.5, if m G (—1, m 2 ) — {mo, mi}, then ( has four complex roots. 

It remains to check the two cases m = mo and m = mi, where the discriminant of 
C vanishes. Since the roots are complex in a neighborhood of each parameter value, one 
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has to check how the roots become repeated as the key parameter value is approached. 
It turns out that in each case, a pair of complex conjugate roots meet on the real axis. 
To see this, we use Descartes' rule of signs to check that the resolvent cubic £ has one 
negative root and a repeated positive root at both m = m and m — mi. It follows that 
one pair of complex conjugate roots has merged into a real root since any other scenario 
would require that zero be a double root of the resolvent cubic. Because the quintic q u (m) 
is squared in the discriminant, the roots return to being two pairs of complex conjugates 
after m passes through these special parameter values. 

To determine whether the repeated real root of ( is positive or negative, we note that 
at either value to = m or m = m±, the shifted quartic rj(w) has the form 

(w — ri) 2 [(w + ri) 2 + r 2 ] = w 4 + {r\ — 2r\) w 2 — 2r\r\ w + r\{r\ + r 2 ). 

Thus, the coefficient of the linear term of rj can be used to determine the sign of the 
repeated real root. Using Sturm's theorem, one can check that this coefficient is negative 
at m — m resulting in a positive value of r±. Then, since r](w) = ((w + c), where c > 
at m — to , we know that the repeated real root of ( when to = to is positive. To find 
this root exactly, we calculate the resultant of ((w) and ^(C) with respect to to. This 
produces a quintic in w given by 

5w 5 - 53w 4 + 98w 3 + 198w 2 + 9w-l. 

The largest root of this quintic, which is approximately 6.9632775, is the positive, repeated 
root of ((w) when to = mo- For the special case to = toi, the derivative of ({w) is a cubic 
with only one real negative root. Consequently, the repeated real root in this case is 
negative. □ 

While it is interesting to discover that the special parameter value to = too has a unique 
solution, it cannot lead to a physical solution to the problem because X\ necessarily has 
to be a different root of ( (since we are excluding the symmetric solutions). But the other 
roots of C are complex. The fact that x\ is complex is also confirmed by solving one of the 
polynomials in the Grobner basis G co i for x\ when to = Too and x 2 = V6.9632775. 

We summarize our findings for the collinear case in the following theorem. 

Theorem 5.2. For the case m e (0,1], there are 12 collinear solutions, one for each 
possible ordering of the vortices. Ifm — 1, all solutions are symmetric and geometrically 
equivalent to the same configuration; otherwise, there are 4 symmetric solutions and 8 
asymmetric solutions. For m e (—1/2,0), there are a total of 6 solutions, 2 symmetric 
and 4 asymmetric, while for m e (—1, —1/2], there are only 2 symmetric solutions and no 
asymmetric solutions. 

Proof. The second polynomial in the Grobner basis G co \, denoted g 2 , is in the elimination 
ideal I4 — I D C[xi,x 2 ,to] and is linear in the variable x±. The coefficient of X\ is found 
to be 8(2to + 1)(to + 2) 2 • q u {m), which is nonzero for to G (—1/2, 1]. By the Extension 
Theorem, we can extend any solution (x 2 ,m) of ( = to a partial solution (xi,x 2 ,m) in 
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V(i4)- Moreover, since g 2 is linear in xi, there is at most one such solution, and X\ must 
be real. 

Next, we check that the value of x\ is distinct from the value of X2 (excluding collision) 
and — X2 (excluding the symmetric solutions.) This is accomplished by substituting x\ = x 2 
into g 2 and then computing the resultant of this polynomial with ((x 2 ,m). A polynomial 
in m is thus obtained and it is easily checked that this polynomial has no roots for m G 
(—1/2,1]. Therefore, the value of X\ obtained by solving g 2 = is distinct from x 2 . A 
similar calculation, using the substitution x\ = —x 2 in g 2 , shows that x\ ^ —x 2 as well. 

At this point, the Extension Theorem can be applied four more times using four of the 
basis polynomials in G co i that are linear in the variables v, u, A and c, respectively, each 
having nonzero leading coefficients. Thus, fixing an m G (—1/2, 1], for each positive root 
of C{ w )i we obtain two possible values of x 2 , each of which extends uniquely to a full, real, 
asymmetric solution of our problem. The precise count on the number of solutions for 



each case then follows directly from Lemma 5.1 and Section 5.1 



The fact that there is one collinear relative equilibrium for each ordering of the vor- 
tices when m > follows from a straight-forward generalization of a well-known result 
in the Newtonian n-body problem due to Moulton [20] • In brief, for each of the n! con- 
nected components of the phase space for the collinear n-vortex problem, there is a unique 
minimum of H restricted to the ellipsoid / = I (see p. 33 of [19] or Section 6.1 of [21] 
for details). Each such minimum is a collinear relative equilibrium. Identifying solutions 
equivalent under a 180° rotation of the plane gives a final count of n!/2. □ 

Theorem 5.3. When m G (—1,0), the signs of the vorticities in a collinear relative 

equilibrium must be arranged as H h (symmetric case only), or as H 1 — or — I h 

(when asymmetric solutions exist). Therefore, when m < 0, it is not possible to have a 
collinear solution where both pairs of vortices with the same strength are adjacent to each 
other. 

Proof. The symmetric case has already been analyzed at the start of this section. For the 
asymmetric case, we note that 

C(l) = -128(1 + 2m) (m + 2) 2 (m + l) 2 

is strictly negative when m G (—1/2,0). We also note that £(0) > and the leading 
coefficient of ((w) is always positive. Since ( has only two positive roots for m G (—1/2, 0), 
it follows that one root is less than one and the other is larger than one. (The two roots 
approach one as m — > — 1/2 + .) The values of x\ and x\ must each be roots of ({w). 
Because we are only considering asymmetric solutions, we have that |xi| < 1 and |a; 2 | > 1 
or vice versa. Therefore, when m G (—1/2,0), the only possible asymmetric orderings of 
vortices have signs arranged asH 1 — or — I h □ 
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6 Asymmetric Strictly Planar Relative Equilibria 



6.1 Eliminating symmetric solutions 

In this section we study strictly planar relative equilibria that do not have a line of sym- 
metry. A major result proved here is that any convex solution with m > or any concave 
solution with m < must contain a line of symmetry. We prove this by saturating the 
Grobner basis in order to eliminate any symmetric solutions, and then showing the result- 
ing system has no real solutions. 

Let J 7 and Q be the Albouy-Chenciner and the unsymmetrized Albouy-Chenciner equa- 
tions in terms of = r?-, respectively, with A' = — 1, Ti = T 2 = 1, and r 3 = T 4 = m ^ — 1. 
Let ecM be the Cayley-Menger determinant written in terms of = rf- and let % be 
the Dziobek equations with A' = —1. The polynomials T^Q^ecM, and % belong to the 
polynomial ring C[m, s 12 , s 13 , s M , s 23 , s 2 4, s u ]- 

We begin by finding a Grobner basis for the ideal I s = (J-,G,ecM,T~i-)- In order to 
accomplish this we first find a Grobner basis Gj g for the ideal J s = (J 7 , Q^cm) an d then 
we compute a Grobner basis Gj s for I s = (Gj g ,T-L). At this stage we saturate with respect 
to the variables S13, su, and S24, to eliminate possible solutions where one of the mutual 
distances has zero length. We also saturate with respect to (S13 — S24), (su— S23), (S13— S14), 



(S23 — S24), (S13 — S23), and (su — S24). Due to Lemma 2.5 and equation (11), saturating 
with respect to these differences is equivalent to eliminating any symmetric solutions. We 
denote the resulting Grobner basis as Gj a . 

Computing Gj s with respect to an elimination order that eliminates all the variables 
except s±2 and S34 yields the following system of two equations in two unknowns: 

s 34 m + S12 - m - 1 = 0, s\ 2 - 2si 2 s 3 4 + s 34 - 1 = 0. 

This system has the two solutions 

, ,/2m + lm\ . , / 1 m + 2\ . . 

(S12, S34 = -7 , S12, S34 = . t 1 — r , 22 

\m + lm+lj \m + 1 m + 1 J 

where one solution is sent to the other one by the transformation m — > — and S12 -H- S34. If 
we compute the Grobner basis of Gj s with respect to an elimination order that eliminates 
all but one Sjj = x, we obtain a polynomial that is the product of the two polynomials 

Pl = 4 (m 2 + 2m + l)x 4 - 4 (5m 2 + 8m + 3)x 3 + 2 (l6m 2 + 21 m + 7): 

- 2 (10 m 2 + 11m + 3)x + Am 2 + 4m + 1, and 

p 2 = 4 (m 2 + 2m + l)x 4 - 4 (3m 2 + 8m + 5)x 3 + 2 (7m 2 + 21 m + 16): 

- 2 (3m 2 + 11m + 10)x + m 2 + 4m + 4, 

where x is one of S13, S14, S23 or S24. For the case when s i2 or S34 is the only variable not 
eliminated, we obtain the products 

(msi2 + S12 - l)(msi2 + S12 - 2m - 1), or 



)x 2 



)x 2 
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(ms 34 + S34 - m)(ms u + s 3i - m - 2), 

respectively. All the Grobner basis computations were performed using Singular 
and Sage [2Tj . 

Analyzing the polynomials p\ and P2 (which we do in the next section), one can prove 
the following: 

Lemma 6.1. For any m G (—1, 1), the polynomial p\ has no real positive roots. For any 
m G (—1, 1), the polynomial P2 has four positive distinct real roots except for m = 0, where 
one of the roots is repeated (at 1). If m G [0, 1], each root of p 2 lies in one of the intervals 
Ji = [0, |] , «/ 2 = [|jl]j^3 = [!)§]' an d ^ 4 = [§>5]- If rn E (—1,0), each root of p 2 lies in 
one of the intervals K x = [0, §], K 2 = [± 1], K 3 = [1, ^±f] and K 4 = [gf , 00]. 



Using Lemma 6.1, we obtain the following fundamental result: 



Theorem 6.2. Let x be a strictly planar relative equilibrium of the four-vortex problem 
with vorticities Ti — T 2 — 1 and T 3 = = m. 

1. For any m > 0, every convex relative equilibrium x has a line of symmetry, and the 
only possible strictly planar asymmetric configurations are concave with the two vortices 
of smaller strength lying on the exterior triangle. Ifm—1, all solutions contain a line 
of symmetry. 

2. For any m < 0, every concave relative equilibrium x has a line of symmetry, and the 
only possible strictly planar asymmetric configurations are convex with equal- strength 
vortices necessarily adjacent. 

Proof. First note that the case A' > is necessarily excluded in the hypotheses of the 
theorem. If m > 0, A' = — A/T < is guaranteed. If m < and A' > 0, then the 
configuration must be convex. Hence, we can assume that A' = — 1. 

If we compute a Grobner basis of Gi s that eliminates all the Sij variables except S12 
and S13, we obtain several polynomials including 

p 3 = (msi2 - 2m + S12 - l)(msi2 + S12 - 1) 

and 

Pi = — 8 ms\ z + 8 (ms 12 + 3 m — s 12 + l)sf 3 — 2 (9 msi 2 — 4 s\ 2 + 14 m — s 12 + 5) s^ 3 
+ 2 (7 msi 2 — 4 s\ 2 + 6 m + Si 2 + 3)si 3 — 3 msi 2 + 2 s\ 2 — 2 m — s\ 2 — 1. 

(23) 

The roots of p 3 are o~\ = 2 ™t^ and a 2 = ^rj. Substituting o"i into p 4 , clearing denomina- 
tors, and rescaling by a constant gives the value pi(s\ 3 ), which is nonzero for m G (—1, 1) 



by the first statement in Lemma [6. 1| Substituting a 2 into p^, clearing denominators, and 
rescaling by a constant gives the value p 2 (si 3 ). For the case m G (—1,1), p 2 has four 
distinct positive real roots. Thus, in order to have a geometrically realizable solution, we 
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must have s\2 = l/(m + 1) and S13 must be one of the four positive roots of p 2 . From 
equation (22), we also have that s 34 = (m + 2)/(m + 1). 
-1, we have l/V- 



2.2 



and 



2.4 



If m £ (0, 1), then 
1 possible differences of 



Since A' = —1, we have l/y/—X' = 1 in Propositions 
S12 < 1 and S34 > 1. Since we have saturated with respect to a 
the remaining four Sij variables, the values of S13, su, S23, S24 must all be distinct roots of 
P2- By Lemma 6.1, two of these values are less than one, and two are greater than one. 



By Proposition 2.4 part 1., the configuration cannot be convex. Since we have saturated 
the Grobner basis to eliminate symmetric solutions, it follows that for m G (0,1), the 
only strictly planar, asymmetric configurations are concave. By Proposition 2.2| part 1., 



vortices 3 and 4 (with equal strength m) lie on the outer triangle, while either vortex 1 or 
2 can lie in the interior of the concave configuration. (The situation is reversed under the 
transformation m h> — .) 

For the special case m = 1, we have that S12 = 3/2 and S34 = 1/2 or vice versa. If 
S12 = 3/2 = <7i, then substitution into p± yields that S13 is a root of p±. But the only real 
roots of pi when m = 1 are 1/2 and 3/2. Thus, we either have s\2 = $13 or S34 = S13 and 



by equation (11), the configuration is necessarily a kite (a symmetric configuration). If 
s i2 = 1/2 = (72, a similar argument with p 2 replacing p\ also yields a kite configuration. 
This completes the proof of part 1. of the theorem. 

If m G (—1, 0), we have S12 > 1, S34 > 1. Taken together with Lemma 6.1, this implies 



that four mutual distances are greater than one and two are less than one. Therefore, from 
Propositions 2J2 and |2.4 it follows that the only real, positive solutions to the system of 
equations given by Gj s correspond to convex asymmetric configurations where the equal- 
strength vortices are adjacent. This proves part 2. of the theorem. □ 



6.2 Proof of Lemma 16.11 



In this section we analyze the polynomials p\ and P2 and prove Lemma 6.1 For pi, we 



use Lemma 3.5 while for p 2 , we make appropriate choices of Mobius transformations. 

Both pi and p 2 are polynomials of degree four. Moreover, if m > 0, their coefficients 
have four sign changes, and thus they have either four, two or zero positive real roots, 
according to Descartes' rule of signs. The roots of the polynomials can be obtained using 
Ferrari's formula for quartic equations [7]. However, we can understand a great deal about 
the solutions by using the resolvent cubic and Mobius transformations. 

One important observation is that if we make the change m 1— > — in the polynomial 
Pi, and clear the denominators, we obtain p^. In particular, this means that the roots of 
P2 for an m G (0, 1) are equivalent to the roots of pi for m G (1, 00). Hence, it suffices to 
study the polynomials for \m\ < 1. 



Discriminants of p\ and P2 

The discriminant of p\ is 256(m+ l) 4 (5m + 3) 2 (2m + l) 2 (m — l) 2 , which is strictly positive 
except for the m- values —1, —3/5, —1/2 and 1 where it vanishes. At m = — 1, p\ reduces 
to a quadratic polynomial with repeated roots at x — 1/2. For m = —3/5, p\ has only 
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complex roots, while for m = —1/2, p\ has a double root at x — 0. At m — 1, there are 
double real roots at a; = 1/2 and x = 3/2. 

The discriminant for p 2 is 256(m + 2) 2 (3m + 5) 2 (m + l) 4 m 2 (m — l) 2 which is strictly 
positive except for the m-values —2, —5/3, —1, and 1 where it vanishes. Focusing on the 
values between —1 and 1, we have that at m = — 1, p2 reduces to a quadratic polynomial 
with repeated roots at x = 1/2. For m = 0, p 2 has a double root at 1 and two other roots 
at (3±V5)/2. At m = 1, p 2 has double real roots at x = 1/2 and x = 3/2. For all values 
of m G (—1, 1), except for m = 0, the four roots of p 2 are distinct. 



Polynomial p\ 



After applying the shift a;4i+(m + l)(5m + 3)/ (4(m + l) 2 ) to remove the cubic term 
in pi, we compute the key coefficients of the shifted quartic to be 



V 



p — 4r 



11m 2 + 6m — 1 

^ — , and 

8(m + l) 2 

(m- l)(5m + 3)(7m + 3) 
16(m + l) 3 



It is straight-forward to check that either p > or p 2 — 4r < (or both) , for each m- value 



in (—1, 1). Using Lemma 3.5, it follows that the roots of p\ are all complex for m G (—1, 1), 
except when m = —1/2, where is a double root. 



We show that 



Polynomial p 2 : The case m G [0, 1] 

Consider the intervals J\ = [0, |], J 2 = [|, 1], J3 = [1, |], and J4 = [|,5]. 
for each m G [0, 1], p 2 has a root in each of these intervals. The result follows from direct 
computation for m = and m = 1, noting that some of the roots are a shared endpoint 
of adjacent intervals. 



rr-interval 



Mobius Transformation 



Ji 
■h 



x 



x 



X 



X 



2 u+1 ' 

W+| 
u+1 ' 

§ u+1 
u+1 ' 

5u+| 
u+1 ' 



m 



m 



m 



m 



a 
a+1 

a 
a+1 

a 
a+1 

a 
a+1 
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Using the Mobius transformations given in the table above, after clearing the denomi- 
nators, we obtain the following four polynomials 

P Jl = - (2 a + l)u 4 - 2 (2 a + l)u 3 + 2 (8 a 2 + 11 a + 4)u 2 + 4 (12 a 2 + 17a + 6)u + 36 a 2 + 48 a + 16 
P, h = 4 a V + 4 (4 a 2 + a) u 3 + 2 (8 a 2 - a - 2)u 2 - 6 (2 a + l)u - 2 a - 1 

Pj., = - (18a + 5)?i 4 - 2(22a + 5)w 3 + 2 (8a 2 - 13a - 2)u 2 + 4a(4a - l)w + 4a 2 
Pj 4 = 4 (3969 a 2 + 3352 a + 704) u 4 + 4 (1764 a 2 + 815 a + 16) w 3 + 2 (392 a 2 - 397 a - 218) u 2 
-2(134a + 45)u- 18a -5. 

Using Descartes' rule of signs, it is straight-forward to show that for any a > 0, each of 
the polynomials above has precisely one positive real root. The values at the endpoints of 
each Ji can be determined by direct substitution. It follows that p 2 = has one solution 
in each of the intervals J±, J 2 , J 3 and J 4 . 

Polynomial p 2 : The case m G (—1,0) 

Consider the intervals K x = [Q,\\,K 2 = [\,l),K z = [1, *±f], and K 4 = [*±f,oo]. We 
show that for each m G (—1,0), p 2 has a root in each of these intervals. 



x-interval Mobius Transformation 



K 1 x=l-^r, m— — \-r 

A 2 x = — rf, m = -r 

z m+1 ' a+1 

X 3 x = "+ 1 , m = — 1- 

° m+1 ' a+1 

K 4 x = u+ m = — 

* m+1 ' a+1 



Pk x = {a 2 + 2a)w 4 + 2(a 2 + 2a)u 3 - 2(4a 2 + 5a + 2)u 2 

- 2(12a 2 + 14a + 4)u - 4(4a 2 + 4a + 1) 
Pk 2 = - 4u 4 + 4(a - 2)w 3 + 2(2a 2 + 7a - 2)w 2 + 6(a 2 + 2a)u + (a 2 + 2a) 
P Ks = 4(4a 4 + 12a 3 + 9a 2 + 2a)u 4 + 8(6a 4 + 17a 3 + 11a 2 + 2a)w 3 

+ 8(2a 4 + 2a 3 - 10a 2 - 9a - 2)w 2 - 8(a 3 + 5a + 2)u - 4a 2 
P Ki = - 16a 3 u 4 - 8(6a 3 + 4a 2 )w 3 - 8(4a 3 + 5a 2 + 2a)u 2 

+ 8(2a 3 + 5a 2 + 2a)w + 4(4a 3 + 12a 2 + 9a + 2). 



Using Descartes' rule of signs, it is straight-forward to show that for any a > 0, each of 
the polynomials above has one positive real root. Again, it is traightforward to check the 
behavior at the endpoints of each interval by direct substitution. It follows that p 2 = 
has one solution in each of the intervals Ki,K 2 , K 3 and K±. □ 
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6.3 The variety of asymmetric configurations 

In this section we study the asymmetric solutions to our problem and show that there are 
exactly eight asymmetric solutions for each m G (—1, 1). We restrict to the case A' = —1 
since A' > only leads to symmetric solutions (as shown in Section 7.3[ ) We first recall 
some definitions and theorems from algebraic geometry (see [6] and [11] for more details). 

Definition 6.3. Let A be a commutative ring. 

1. An ideal I of A is said to be real if, for every sequence a±, . . . , a p of elements of A, we 
have 

a\ + . . . at G I ==>■ di G I, for i = 1, . . . , p. 

2. \/J is the smallest real ideal of A containing / and is called the real radical of the 
ideal /. 

Let A; be a field. If I C k[x\, . . . , x n ] is an ideal, we denote by V(I) the set 

V(J) = {(ai, . . . , a n ) G k n : f(a u . . . , a n ) = for all / G /}. 

V(J) is an affine variety. In particular if / =< fx, . . . , f s >, then V(J) = V(/i, . . . , f s ). 
If k is the field of the real numbers M we will say V(7) is a real algebraic variety. Note 
that while this terminology is common in algebraic geometry books, it is different from 
the terminology frequently used in real algebraic geometry (see [6], for example). 

Definition 6.4. Let V C k n be an affine variety. Then we set 

I<y) = {/ e fcfci) ■■■■, x n) ■ /(ai, . . . a n ) = for all (ax,... a n ) G V}. 
1(V) is an ideal and it is called the ideal of the variety V. 

We are now ready to state a version of the Real Nullstellensatz: 

Theorem 6.5 (Real Nullstellensatz). Let k be a real closed field and I an ideal of k[x\, . . . , x r , 
Then I(V(/)) = ^fl . 

The dimension of an affine variety V C k n , denoted dimK, is the degree of the affine 
Hilbert polynomial of the corresponding ideal 1(V) C k[xi, . . . [TT] for more 

details). This degree can be easily computed using Singular. In the case k = M. one 
needs to know the ideal of the variety, or by the Real Nullstellensatz, the real radical of the 
ideal. The real radical can be computed using the realrad. lib [26] library of Singular 
[T2] . More details about the algorithms can be found in the paper [25] . We now apply the 
above theory to our problem. 

Theorem 6.6. Consider the ring M[si2, si3> su, $23) S24, S34, m]. The asymmetric relative 
equilibria configurations form a one- dimensional real variety V C M7 . 
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Proof. Consider the Grobner basis Gj s for the ideal I s and saturate with respect to 

Sl3; S14, S24, (S13 — S24), (S14 — S23), (S13 — S14), (^23 — S24), (^13 ~ s 23)) an d (^14 ~ ^24), 

as before. In addition, saturate with respect to (msx2 — 2m + S12 — 1) to eliminate the 
case s±2 = (2m + l)/(m + 1) for which the solutions are complex. We obtain the following 
polynomials: 

fl — S13 + S14 + S23 + S24 — 2S34 — 1 

h = s u - s 34 + 1 

fa = s 3im + s 34 - m - 2 

/4 = 2S24 — S14S34 — S23S34 — 4S24-534 + 2S3 4 + 2s 2 3 
U = 2^23^24 - 2S 23 - 2S 2 4 + S34 

fe — 2S14S24 — S34 

f 7 = 2sl 3 + S14S34 - 3S23S34 + 2s 2 4 - S34 

f S = 2si 4 S23 - S14S34 - S23S34 + S34 

2 

14 



fg — 2s\ 4 — 3S14S34 + S23S34 — 2^14 — 2s 2 3 — 2^24 + 3^34 + 2. 



Let / =< fx, . . . , fg >. We want to find the dimension of the real variety V(J). First, 
we find I(V(7)) or yfl. This can be computed using the realrad.lib [2S] library of 
Singular [T2]. However in this case it turns out that \[l = I, and hence fx, . . . , fg are 
the generators of the real radical of /. The degree of the affine Hilbert polynomial of 
^1 = 1 is one, and thus V(7) is a one dimensional real algebraic variety. □ 

We now want to see if the variety contains singular points. First we recall some 
definitions and theorems. 

Definition 6.7. Let V C k n be an affine algebraic variety then the Zariski tangent space 
of V at p — (px, ■ ■ ■ ,p n ) denoted by T p (V) is the variety 

T p (V)=V(d p (f):feI(V)) 

where d p (f) = ^=1 §^(p)(^i ~ Pi)- 

Clearly if I(V) =< fx, . . . J s > then T P (V) = V(d„(/i), . . . , d p (f s )), and it is the 
translate of a linear subspace of k n . Recall that if V C k n is an affine variety, then V is 
irreducible if and only if I(V) is a prime ideal. Thus we have the following definition: 

Definition 6.8. Let V C k n be an irreducible variety. A point p in V is nonsingular 
(or smooth) provided that dimT p (V) = dim V. In other words, if I(V) =< fx, ■ ■ ■ , f s >, 



then p is nonsingular if and only if the rank of the Jacobian matrix J 
n — dim(V). Otherwise, p is a singular point of V. 



dfi 



is equal to 



The notion of Zariski tangent space is defined even at a singular point. However, when 

V C W 1 is an irreducible variety and p e V is a nonsingular point, a neighbourhood of p in 

V is a l7°° submanifold of 1R™ and the tangent space of V at z (in the l7°° sense) coincide 
with the Zariski tangent space. If p is singular, then the dimension of T P (V) is bigger than 
the dimension of V. 
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Theorem 6.9. The variety of the asymmetric strictly planar relative equilibria configura- 
tions has no singular points and hence it is a smooth one- dimensional manifold. 



Proof. The Jacobian matrix of (/i, . . . , fg) is 



"0 


1 


1 


1 


1 


-2 





1 














-1 




















m + 1 


S34 - 1 








-S34 


-S34 + 2 


4 («24 - S34) 


-814 ~ S23 - 4(fi24 - S34) 














2 s 24 - 2 


2s 23 -2 


1 











2s 2 4 





2S14 


-1 











S34 


4 S23 - 3 s 34 


2 


Sl4 - 3s 2 3 — 1 











2 S23 - S34 


2 S14 - S34 





S14 ~ S 2 3 + 1 





.0 





4 S14 — 3 S34 — 2 


s 34 - 2 


-2 


—3 S14 + S23 + 3 






Using Gaussian elimination we obtain the following matrix 



where 



1 












-1 


" 


1 1 


1 




1 




-2 





s 34 


4 S23 - 3 s 34 


2 


S U - 


3S23 - 


1 





2s 2 4 - 


2 


2S 23 - 


2 


1 































A 






















































Qi 


Q2 
















m + 1 


S34 — 2 








A = 





a 
















b 
















c 










and a, b, c are polynomials that belong to the ideal of the variety (and hence vanish on the 
variety), and 

qi = -4(s 23 - s 2 4)(2s 23 + 2s 24 - 2s 34 - 1), 

q 2 = -8 s 23 s 2 4 ~ 8 s 24 2 + 8 S24 S34 + 6 s 24 + 4 s 23 - 4 s 34 . 

Since s 34 7^ 0, s 2 4 — 1 7^ 0, q\ 7^ 0, and m 3 — 1 and s 34 — 2 are not both zero on the variety, 
it follows that the rank of the Jacobian matrix is 6. Since dim(V r ) = 1, by Definition 6.8 
the variety has no singular points. It follows that it is a smooth manifold. □ 



Theorem 6.10. For each value of m in (—1,1), there are exactly four (eight counting 
reflected solutions ) strictly planar asymmetric relative equilibria configurations. They are 
convex if m < and concave if m > 0. 
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Figure 2: Projection of part of the manifold of asymmetric configurations onto the plane 
(m, S23). As m — >■ —1, S23 — > 00 along two branches while S23 — > \ along the other two 
branches. 



Proof. The proof consists of showing that the system of equations f\ = 0, . . . , fg = has 
four solutions for each value of m in (—1,1). The convexity of the configuration then 



follows from Theorem We use an algorithm for the triangular decomposition of semi- 
algebraic systems. Such an algorithm, given a system of equations and inequalities S, 
computes simpler systems S%,...,Sk such that a point is a solution of the original system 
5* if and only if it is a solution of one of the systems Si, ...,Sk- Each of these systems has 
a triangular shape and remarkable properties: for this reason it is called a regular semi- 
algebraic system and the set of the S%, ...,Sk is called a full triangular decomposition of S. 
See [U] and references therein for some background. The algorithm is detailed in [U] and 
it is available in Maplel5™ via the RealTriangularize command of the RegularChains 
package. A full triangular decomposition of the system f\ = 0, . . . , fg = is given below. 
If m = we have the following two systems: 



sia - 1 = 
S13 + s 23 - 3 = 
1 = 



S14 

So, - 3s 2 3 + 1 = 
1 = 



'23 

S34 ~ 2 = 



Sia - 1 = 

s 13 - 1 = 

S14 + S24 — 3 = 

S23 ~ 1 = 



'34 



3^24 







S34 - 2, 
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if m — 1 we have: 



2s 12 -1 = 


2s 12 -1 = 


2s 13 -1 = 


2si 3 -3 = 


2s M -1 = 


2s 14 -3 = 


2s 23 -3 = 


2s 23 -1 = 


2s 24 -3 = 


2s 24 -1 = 


2s u -3 = 


2s 34 -3 = 



and if m G (—1, 0) U (0, 1) we obtain: 



su - s 34 + 1 = 
Sl 3 + s 14 + s 23 + s 24 
s 34 Si4 + (s 34 — 2)s 23 — ^c 24 
(2s 24 - 2)s 23 - 2s 24 + s 34 = 

Ast + (-8s 34 - 4)s\ A + (4s 2 4 + 6s 34 + 4)s 2 24 + (-4s^ 4 - 2s 34 )s 24 + 



2s 34 -1 = 
+ 4s 24 s 34 



2s 2 
Z5 34 



"24 

(m + l)s 34 







This decomposition immediately shows that the system has no real solutions if m = — 1, 
two real solutions if m = 1, and four real solutions if m = 0. It remains to study the 
case m e (—1,0) U (0, 1). Solving the last equation of the last system of the triangular 

Substituting this into the preceding equation of the same 



m+l ' 



decomposition yields s 34 

system, after clearing the denominators, yields p 2 (s 2 4). From Lemma 6.1 we know that p 2 
has four distinct positive real roots for m G (—1,0) U (0, 1). If we fix one such solution 
for s 24 and substitute into the first four equations of the system, the resulting system is 
Ax = b, where x = (s 12 , s 13 , s u , s 23 ), 



A 



10 
11 1 

s 34 (ss4-2) 

2(s 24 -l) 



and b = [—1 + s 34 , — s 24 + 2s 34 + 1, s 24 — 4s 24 s 34 + 2sg 4 , 2s 24 — s 34 ] T . Since the matrix A 
is invertible, Ax = b has exactly one solution for every vector b. This proves the theorem 
since it implies that the system fi = 0, . . . fg = has four real positive solutions. 

□ 



Remark. Note that for m = 1 we obtain two solutions which correspond to an equilateral 
triangle configuration with a vortex at the center. These are symmetric configurations that 
have not been eliminated because in the case of equal vorticities, there is an additional 
symmetry that we did not exclude when saturating. Also, note that the polynomials 
/i, . . . , fg were obtained saturating with respect to (msi 2 — 2m + si 2 — 1). This eliminated 
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the real solutions of the system f x — 0, . . . , fg — for m 6 (— oo, —1) U (1, oo). However, 
we know that the number of real (positive solutions) of the original system (e.g., before 
saturating) for m G (— oo, — 1) U (l,oo) is the same as the number of solutions of fx = 
0,...,/ 9 = 0forme (-1,1). 



7 Symmetric Strictly Planar Relative Equilibria 

In this section we investigate all possible symmetric, strictly planar, relative equilibria 
in the four-vortex problem with Ti = I^ = 1 and T3 = r 4 = m. The four possible 
configurations are an isosceles trapezoid, a concave kite, a convex kite, or a rhombus. We 
separate the kite configurations into two cases based on the sign of A'. 



7.1 The isosceles trapezoid family 



Without loss of generality, suppose that the vortices are ordered sequentially around an 
isosceles trapezoid, so that the lengths of the diagonals are given by r\% = r 24 and the 
congruent legs have length ru = T23- Due to the symmetry of the configuration, the four 
vortices lie on a common circle, that is, the isosceles trapezoid is a cyclic quadrilateral. 
By Ptolemy's theorem, we have that 



13 



24 



^12^34 + r 14 . 



(24) 



^24, "^14 = "^23 and equation (24) will also 



A choice of six mutual distances satisfying r i3 
satisfy ecu = 0. 

In contrast to the Newtonian four-body problem, the isosceles trapezoid family of 
four-vortex relative equilibria can be solved completely (by hand) in terms of the vortex 
strength m (compare with [IH])- Let x = r^/r \i and y = r^/r 12- Rather than assuming 
a specific scaling on the distances (e.g., A' = — 1), we will solve for x and y in terms of m. 



Using the symmetries of the configuration and equation (24), equation (10) reduces to 

2(y 2 - x 2 )(y 2 - 1 



x(2y 2 



x 



1) 



0. 



(25) 



This relation between x and y is both necessary and sufficient for the trapezoid to be a 
relative equilibrium. 

Focusing on the first and last columns of equations (|9]), the first and last equations are 
easily satisfied due to symmetry. The middle four equations in pi) are equivalent given 



that (25) holds. They determine a formula for m given by 



m 



x(l-y 2 ) 



X' 



(26) 



where we have used | 
of radius r c 
equation 



r jk'Tjirki/{4:r c ) for the area of a triangle circumscribed in a circle 
Solving equation (26) for y 2 and substituting into equation (25) leads to the 



x(l - x)(l + x)(x 2 (2m + 1) - m(m + 2)) 



0. 



(27) 
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Since x represents the ratio of two distances, there are only two possibilities. If x = 1, 



then equation (25) quickly gives y — 1 and the configuration is a square. The second 



equation in <\9h then gives that m = 1 and all vortices must have the same strength. The 



other possibility is that x 2 = m(m + 2) / (2m+ 1). Substituting this value into equation (26 ) 
and solving for y 2 yields 




, m(m + 2) 
m + 2 - \ — ± — —-t 
1 2m + 1 

Note that the necessary condition x 2 > is satisfied only for m > and —2 < m < —1/2. 
In order for the expression for y 2 to be real, we must be in one of these two regimes. 
However, while y 2 > holds for m > 0, it is not satisfied on —2 < m < —1/2. Therefore, 
there does not exist an isosceles trapezoid solution for the case m < 0. 

Theorem 7.1. There exists a one-parameter family of isosceles trapezoid relative equilibria 
with vortex strengths T\ — Y 2 — 1 and r 3 = T 4 = m. The vortices 1 and 2 lie on one base 
of the trapezoid, while 3 and 4 He on the other. Let a = m(m + 2)/ (2m + 1) . If r 13 = r 2 4 
are the lengths of the two congruent diagonals, then the mutual distances are described by 

\ 2 / \ 2 i 

— I = a 

ri2 



(S =5(" + *-V5) 

and ( — J = ^(m + 2 + A/a). 
\ r i2 / 2 



T/izs family exists if and only if m > 0. T/ie case m = 1 reduces to the square. For m/1, 
i/ie larger pair of vortices lie on the longest base. 



Proof. The formula for r 2 3 /r 2 2 comes from equation (24). Substituting m = 1 into the 
formulas above quickly gives r±2 = r i4 = r 34 = r 23 and r±3 = r 24 = \/2ri2, a square. The 
only other item that remains to be shown is the fact about the larger vortices lying on the 
longer base. This follows since a > 1 if and only if m > 1. □ 

Remark. 1. If < m < 1, r±2 is the longest base length and the formulas in the theorem 
give r 34 < ri4 < ri 2 < ri 3 . On the other hand, if m > 1, then r 34 is the longest base 
length and we deduce ri 2 < r i4 < r 34 < r i3 . Both cases agree with the conclusion of 



statement 1 in Proposition 2.4 



2. The case m > 1 is identical to the case < m < 1 through a rescaling of the vortex 
strengths and interchanging bodies 1 and 3, and bodies 2 and 4. Specifically, replacing 
m by 1/m and interchanging distances r 12 and r 34 leaves the three equations for the 
ratios of mutual distances unchanged. The vortex strengths map to Ti = T 2 = 1 and 
r 3 = r 4 = l/m<l under this transformation. 

3. As m —?• 0, a — >■ and consequently r 34 — > 0. The limiting configuration is an equilateral 
triangle with vortices 3 and 4 colliding. As m increases, so does a, and the smaller base 
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length r 34 approaches the larger one r 12 as m — > 1. The ratio of the diagonal to the larger 
base also increases monotonically with m. However, the ratio of the leg to the larger 
base begins at 1 {m — 0) and decreases to a minimum value of approximately 0.904781 
before increasing back to 1 at the square configuration. This minimum value occurs at 
m ~ 0.234658, which is the only positive root of the quartic 8m 4 + 19m 3 + 9m 2 + m — 1. 
Just as with the Newtonian case (see [ID]), the range of y is surprisingly small, confined 
to the interval [0.904781,1]. 



7.2 Kite configurations: A' < 

In this section, we consider kite central configurations when A' < 0. Such configurations 
contain two vortices that are symmetrically located with respect to an axis of symmetry 
while the remaining two vortices lie on the axis of symmetry. The configuration formed by 
the vortices can either be concave or convex. In this section, we focus on configurations 
with only one line of symmetry. Those configurations with two lines of symmetry, i.e., the 



rhombus configurations, are studied in detail in Section \7A 

Our goal is to count the number of kite configurations as m varies and describe the type 
of possible configurations. Recall that J 7 , Q, cqm and % represent the Albouy-Chenciner 
equations, the unsymmetrized Albouy-Chenciner equations, the Cayley-Menger determi- 
nant and the Dziobek equations, respectively, in terms of the variables Sij = rfj, with 

A' = —1. Denote the complete system of these equations as E. Our goal is to prove the 
following theorem. 

Theorem 7.2. For each value of m G (— oo, — 2) U (— 1,0) U (0, oo), system £ has ex- 
actly four (real, positive) solutions that correspond to symmetric kite configurations. If 
m G (0, oo), such configurations are concave. If m G (— 1,0), then there are two convex 
configurations with T 3 and T 4 on the axis of symmetry, and two concave ones with Ti and 
r 2 on the axis of symmetry. If m G (— oo, —2), the situation is reversed and there are two 
concave configurations with T 3 and T 4 on the axis of symmetry, and two convex ones with 
Ti and T 2 on the axis of symmetry. There are no other strictly kite configurations with 
A' < for other values of m. 

Remark. 1. As we show below, the case m — 1 is considerably different from the New- 
tonian four-body problem. In fact, it is known for the Newtonian case that there are 
two distinct symmetric kite configurations for four equal masses: an equilateral tri- 
angle with a mass in the center and an isosceles triangle with a mass on the axis of 
symmetry [TJ. Instead, for four vortices of equal strength, the only concave kite central 
configuration is an equilateral triangle with a vortex at the center of the triangle. 



2. Based on this theorem and Theorem 6.10 we conclude that there are no concave relative 
equilibria when m G (—1,-1/2). This is another contrast with the Newtonian four- 
body problem, where it is shown in [H] that for any choice of four positive masses, 
there exits a concave central configuration. 
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3. The exact count on the number of concave symmetric solutions given in Table [T] is 
precisely twice the values expressed in the theorem. This is due to the fact that each 
solution found can be reflected about the axis of symmetry, a transformation that does 
not change the mutual distances r^, but does alter the positions 

We have two possibilities, either vortices 3 and 4 lie on the axis of symmetry or vortices 
1 and 2 do. It is convenient to study the two cases separately, since the conditions imposed 



by the symmetry are different. The first case will be explored in Lemma 7.3 , and the second 



in Lemma 7_A The proof of Theorem 7/2 follows immediately from the two lemmas. 

Let us consider the first case. If vortices 3 and 4 lie on the axis of symmetry, then 
•S23 = S13 and S24 = su. We compute a Grobner basis for the ideal (F, Q, ecM, H), where 

G, &cm, 'H are obtained from J 7 , Q, com and H, respectively, by imposing the conditions 
S23 = Sis and s 2i = sy. Then we saturate with respect to Sj 3 — su, in order to exclude 
the rhombus configurations. We also saturate with respect to si 2 , S13, su, S34 and m. This 
computation yields the following ten polynomials: 

9\ =S3<im + S13 + s u + S34 - m - 2 

92 =Si^m + s 14 m + s 12 + s 13 + s u + s 34 - 3m - 3 

gz =2si 2 s 34 - 2si 3 s 34 - 2si 4 s 34 + 2s 2 34 - s i2 + 2si 3 + 2s 14 + s 34 - 3 

9i =2Si 3 Si4 - S13 - S 14 + S34 

95 = 2 S 13 + 2s^ 4 — 2s 13 S34 — 2^14534 + S12 — 4s 13 — 4Si4 + S34 + 3 

g 6 =2si 2 si 3 + 2si 2 si 4 - si 2 m - 2s 12 - s 13 - s u + s 34 

g 7 =si 2 m 2 + 2sl 2 + 2si 3 s 34 + 2si 4 s 34 - 2sg 4 - 3si 2 m - As 12 

- Si3 - su - 2s 34 + 2m + 4 
g 8 =2s\ A m + 2si 2 si 4 + 2s^ 4 + 2si 4 s 34 - Qs u m - si 2 + si 3 - 5si 4 - s 34 

+ 2m + 1 

99 = 4s 44 — 4s^ 4 S34 + 2Si 2 Si4 — 10s 44 — 2S13S34 + 2S14S34 + 2Sg 4 — S12 



+ s i3 + 7su + 2s 34 - 3 

I s i2 s i4 — 2si2Si4m — 4si2Si4 — 2s^ 4 — 2S13S34 + 2s^ 4 + Si 2 m + S12 
+ 2si3 "I" 2si4 + S34 — 3. 



Let us denote the corresponding system of equations gi = 0, . . . , gio = as G. Then the 
positive solutions of G give the symmetric kite configurations described above. An exact 
count of these solutions is described in the following lemma. 

Lemma 7.3. For each value ofm G (0, 1), system G has exactly four real positive solutions 
corresponding to concave configurations. For each value of m G (—00, —2) U (— 1,0), 
system G has exactly two real positive solutions. Such solutions correspond to convex 
configurations if m G (— 1,0) and to concave configurations if m G (—00, —2). Form = 1, 
system G has exactly two real positive solutions corresponding to an equilateral triangle 
with either vortex 3 or 4 at the center. There are no positive real solutions for other values 
of m. 
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Proof. Computing a full triangular decomposition of system G with the positivity condi- 
tions S12 > 0, Sis > 0, S14 > and s 34 > 0, we obtain five simpler systems. If m — 1, we 
have the following two systems: 



If m = 0, we have the following two systems: 



2S12 


-3 


= 


2S13 


-3 


= 


2s M 


- 1 


= 


2s 34 


- 1 


= 0, 




(28) 



sia - 1 = 

Si 3 - 1 = 

4s 14 -1 = 
4s 34 -3 = 0, 



T 4 




lime (-oo, -2) U (-1/2, 0) U (0, 1), we obtain the system 



n = 



512 + (— m 2 — 2m)s 34 — 1 + m 2 = 

513 + S14 + (m + l)s 3 4 - 2 - m = 
2s 2 4 + ((2m + 2)s 34 - 4 - 2m)si 4 + 



-2 - m)s 34 + m + 2 = 



(4 + 6m + 2m 2 )sg 4 + (—3 — 3m 2 — 6m)s 34 + 2m + m 2 



0. 



For all other m-values the simpler systems have no real solutions. 

The case m = 1 is straightforward, since T\ and T2 each have a unique solution corre- 
sponding to an equilateral triangle with either vortex 3 or 4 at the center. The case m = 
is uninteresting and we do not analyze it. It remains to study the solutions of system T 5 . 

• Solutions of T 5 with S34 > and s u > 

The first two equations of T5 are linear in s\2 and si 3 , so that to each solution (si 4 , s 34 ) 
of the last two equations correspond one solution (S12, S13, S14, S34) of T5, provided that 
the determinant of the matrix of the coefficients of the linear system is not zero. Hence, 
we focus our attention on the last two equations of T 5 . Let A = 2(2 + 3m + m 2 ), B = 
— 3(m + l) 2 , and C — m(2 + m). Then the last equation of T 5 can be written as 



5i(ss4) 



As 2 4 + Bs u 



C = 



and can be viewed as a parametric equation of a parabola. Clearly B < for all m 7^ — 1, 
and A > for m E (—00, —2) U (—1, 00). Thus A > 0, B < on the domain where T 5 is 
valid. The vertex of the parabola has coordinates 



2^' gi ("2^) 



3(m + 1) 



1 (m - l)(m 2 - 4m - 9) 



2(2 + 3m + m 2 )' 8 



m + 2 
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It is easy to see that — ^ > and qi(-^j) < in the interval of interest, and hence the 
equation has two real solutions for each value of me (—00, —2) U (— 1/2, 0) U (0, 1). Since 
C > for m G (—00, —2) U (0, 00) and C < for m € (—2, 0), it follows that the equation 
has two positive solutions for m G (— 00, —2) U (0, 1) and only one positive solution for 
m G (-1/2,0). 

Let us denote as /3i and (3 2 (with 0i < f3 2 ) the (real) solutions of (71(534) = 0. Then 
through a standard analysis it is possible to find bounds for f3 1 and (3 2 . Such bounds are 
summarized in the table below. 



solutions 


m G (—00, —2) 


mG (-|,0) 


m G (0,1) 




< ft < \ 


ft <o 


< ft < \ 


P2 


&>1 


< ft < 1 


< ft < 1 



Now consider the third equation of system T5. Let A' = 2, B' = (2(m+ l)s34 — 4 — 2m), 
and C" = (m + 2)(1 — s 34 ). Then this equation takes the form 

q 2 (s u ) = A's 2 u + B's u + C' = 0. 

Since this is a quadratic equation in su, there are going to be two solutions for each s 34 
and m. We will denote the (real) solutions corresponding to S34 = (3\ as (3\ and (3j (with 
Pi < Pi), an d the ones corresponding to s 34 = f3 2 as ^ and /3 2 (with ^ < 13%). The vertex 
of this parabola has coordinates 



B' 



2A 



-,,<12 



B' 



2A' 

(m + l)(l -S34) + 1 1 



m 



9 v m + 1) 2 4 + m(m + 2)s 34 - y (m + 2) ) . (29) 



We can then view the ordinate of the vertex, denoted as r(s 34 ), as a parabola with coeffi- 
cients in m: 

r(s 34 ) = --(m + 1) 2 4 + m(m + 2)s 3 4 --{m + 2)= A"s 2 34 + B"s u + C" = 0. 



The coordinates of the vertex of r(«34) are 



B" 
'2A' [ 



B 1 ' 



2A" 



m{m + 2) 1 m(m + 2) 
(m + l) 2 ,_ 2 (m + l) 2 



CASE I: m G (0, 1) 

Since > 0, r (-^7) < 0, and A" < for m G (0, 1), it follows that q 2 (-^7) < 

for m G (0, 1). Consequently, since C > in the same interval, g 2 = has two positive 
solutions corresponding to each solution of q± — for m G (0, 1). In other words, there 
are four (positive) possible values for su, namely (5\ > 0, /3f > 0, (5\ > 0, and /3 2 > 0. 
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CASE II: m G (-§,0) 

For m G (— 1,0), we have A' > and C > 0. This implies that corresponding to 
s 34 — P2, there are either two real positive solutions or no real solutions. We want to 
show that there are two positive real solutions. Clearly —^77 < 0, r (—jX") > 0, A" < 
and C" > for m G (— \, 0), so that r = has one positive and one negative root. If we 
show that the positive root is less than s 34 = /3 2 , then it follows that q 2 (—^7) < 0, and 
consequently q 2 = has two positive solutions when s 34 = f3 2 (i.e., fi\ > and f3 2 > 0). 

Subtracting f3 2 from the positive root of r = yields 

— (— m 3 — 7 m 2 — 4im + ym — 7 m — 8 a; + y + 3) 
4(m + 2) (m+ l) 2 

where x = a/— m(m + 2), and y = ^J(m — l)(m + l)(m 2 — 4m — 9). Numerically it is 
easy to see that A < for m G (— |, 0). However, we proceed more rigorously. Taking the 
numerator of the expression above and squaring the expressions for x and y yields three 
polynomials. Computing a lex Grobner basis for these three polynomials, we obtain 

m(5m 3 + 7m 2 + 3m + 9)(m + 2) 2 (m + l) 2 

as one of the basis elements. It is trivial to see that this polynomial has no roots in [— |, 0], 
that A(0) < 0, and that A(— |) < 0. Since A is continuous on [-1/2,0], we deduce that 
A < for m G (— |, 0), as desired. 

CASE III: m G (-00, -2) 

If s 34 = P 1 and m G (-00, -2), then C < 0. Since A 1 > 0, it follows that q 2 = 
has one positive and one negative real solutions corresponding to If S34 = f3 2 and 
m G (-00, -2), then C > 0. Since r (-7JJ) < and A" < it follows that q 2 (-^7) < 0. 
This together with the fact that — ^7 > 0, A' > and C > implies that, corresponding 
to (3 2 , q 2 = has two real positive solutions, (i.e., $\ > 0, and $\ > 0). 

• Solutions of T 5 w?£/i s 34 > ; s u > 0, and s 13 > 

We now study which solutions among the ones found above have Si 3 > 0. The resultant 
of the second and third equations of T 5 with respect to s 12 is 

R = 2s 2 3 + ((2m + 2)s 34 - 2m - 4)s 13 + (-2 - m)s 34 + m + 2. 

If si 3 is replaced with si 4 , this polynomial transforms into the one that appears in the 
third equation of T 5 . Hence our discussion regarding s i4 applies to s i3 . Note that since we 
saturated with respect to Si 3 — Si 4 , every solution must have Si 3 7^ Si 4 . Since i? = has 
two solutions for each m and S3 4 , when si 4 is one such solution, S13 must be the other. 

• Solutions of T 5 wz^/i s 34 > 0, Si 4 > 0, S13 > 0, and s 12 > 

We now study which solutions among the ones found above have Si 2 > 0. Consider 
the first equation of T 5 and rewrite it in the form 

s 12 = (m 2 + 2m)s 34 + l -m 2 . (30) 
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s 34 




Figure 3: This figure summarizes the analysis of the sign of si 2 . The grey areas represent 
regions where s\2 > 0. The white areas represent regions where Si 2 < 0. Along the solid 
lines S12 = 0. 

This equation has a unique solution for each value of S34 and m. Solutions corresponding 
to S34 = Pi will be denoted aij. We can view this expression as defining a function of the 
two variables m and S34. The sign of this function, i.e., the sign of Si 2 , is summarized in 
Figure [3j As mentioned earlier, the last equation of T 5 gives that s 34 = (3i or s 34 = /3 2 , 
and when /3j > 0, we have < S34 < 1 for m G (— 1,0) U (0, 1). Hence Si 2 > for any 
positive solution of qi(su) = 0. 

The situation is a bit more delicate when m G (—00,— 2). In this case, 02 > 1 ; an d 
thus we cannot arrive at a definitive conclusion just looking at Figure [3j However, using 
standard methods, one can prove that the following function of m, defined by 

d(m) = f3 2 (m) y — — s , 

is always positive (where the graph of /(to) = ^^"+2) ^ s ^ ne boundary of the grey region 
in Figure [3] in the interval (—00,— 2)). Thus, (3 2 lies in the region where s\2 > when 
m G (—00, —2). Consequently, there are always positive solutions of T 5 corresponding to 
/3 2 . On the other hand, < (5\ < \ for m G (—00, —2), and thus there are no positive 
solutions of T 5 corresponding to (5\ in this interval. 

In summary, we have shown that the solutions (si 2 , S13, S14, s 34 ) T with all components 
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positive have the form 











a 2 




a 2 


Pi 




PI 




PI 




PI 


PI 


1 


PI 


1 


PI 




PI 


Pi 




Pi 




P2 




P2 



for m G (0, 1), while they have the form 



a 2 




a 2 


PI 




PI 


PI 


•> 


PI 


A 




A 



for m G 



-00, -2) U (-|,0). 

We are left with deciding which solutions are convex and which are concave. For 
m G (0,oo), it follows from the main result in [3] that, if the kite is convex, it must be 
a rhombus. Since the rhombus case was excluded, these solutions correspond to concave 
configurations. For m G (— 1,0), we have that s 34 = (3 2 < 1. However, according to 



Proposition 2.2, since A' = —1, the interior side connecting equal vorticities in a concave 
configuration is greater than 1. Hence, the configurations are, in this case, convex. For 
m G (—00, —2) we have that S34 = (3 2 > 1. Furthermore, substituting s 34 = f3 2 {m) 



into equation (30), we obtain an expression that is function of m alone. Examining this 
function, it is easy to show that si 2 = a 2 {m) < 1 for all m G (— 00, 0). Thus, from 



Propositions |2.2| and |2.4[ it follows that the solutions correspond, in this case, to concave 
configurations. 



This completes the proof of Lemma 7.3 



□ 

We now consider the other case, where vortices 1 and 2 lie on the axis of symmetry. We 
compute a Grobner basis for the ideal , Q', e' CM , H'), where T' y Q', e' CM , %' are obtained 
from T,G,ecM an d H, respectively, by imposing the conditions s 14 = Si 3 and s 24 = s 2 3- 
We saturate with respect to Si 3 — s 23 in order to exclude the rhombus configurations. We 
also saturate with respect to si 2 , S13, S23, S34 and m. This computation yields a polynomial 
system that we will call G%. Then we can prove the following lemma. 

Lemma 7.4. For each value of m G (l,oo), system G\ has exactly four real positive 
solutions corresponding to concave configurations. For each value of m G (— 00, —2) U 
(— §,0) ; system G\ has exactly two real positive solutions. Such solutions correspond to 
concave configurations if m G (— 1,0), and to convex configurations if m G (—00, —2). For 
m = 1, system G\ has exactly two real positive solutions corresponding to an equilateral 
triangle with either vortex 1 or 2 at the center. There are no positive real solutions for 
other values of m. 
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Proof. Computing a full triangular decomposition of Gi, with the positivity conditions 
S12 > 0, S13 > 0, Sx4 > and s 34 > 0, we decompose G\ into three simpler systems Si, S 2 
and ^3. For m = 1, we have two systems, Si and S 2 . They are linear, and each of them 
has a unique solution corresponding to an equilateral triangle configuration with vortex 1 
or 2 at the center. For m G (—00, —2) U (— |, 0) U (1, 00), we obtain a third system S3. This 
last system, possibly after computing a lexicographic Grobner basis, can be reduced to T 5 
by performing the transformation (m, Si 2 , s i3 , s 23 , s 34 ) (^, s 34 , Si 4 , Si 3 , Si 2 ) and clearing 
the denominators. This transforms each of the four polynomial equations in S3 to the 



equations of T 5 . Consequently, these cases follow from Lemma 7.3 For all other m- values 



the simpler systems have no real solutions. □ 



Remark. 1. The proofs of Lemmas 7.3 and 7.4 yield some interesting information on the 



type of kite central configurations. For example, when m G (0, 1), there are solutions 
only when the smaller strength vortices (3 and 4) are on the line of symmetry. In this 
case, there are two geometrically distinct configurations for each choice of the central 
vortex. 

2. As described in the introduction, the equilateral triangle with a central vortex bifurcates 
into both concave kite configurations and asymmetric concave configurations as m 
decreases through m = 1. The fact that this particular solution shows up in both 



Sections and is hardly surprising. In fact, when solving the Albouy-Chenciner 
equations for m = 1 without any restrictions on the variables, this solution occurs 
four times with a multiplicity of four. The Hessian matrix D 2 (H + A J) evaluated at 
this solution has a null space of dimension 2 (excluding the "trivial" eigenvector in 
the direction of rotation). One can check that there are eigenvectors in the null space 
corresponding to both a symmetric and asymmetric perturbation of the solution. In 
fact, the symmetric perturbation gives rise to the concave kite configurations described 
above. The central configuration consisting of four equal vortices on an equilateral 
triangle with a vortex at the center is a highly degenerate solution. 

7.3 Kite configurations: A' > 



In Section [2TT] we observed that, if the vorticities have different signs, we could have central 
configurations with A' > 0. In this case, taking Ti = T 2 = 1 and T 3 = T 4 = m < 0, the 
configuration is necessarily convex with vortices 1 and 2 on one diagonal and vortices 3 
and 4 on the other. In this section, we first show that these convex solutions must contain 
a line of symmetry. Based on the main result in [3] (applicable to our problem only 
when m > 0), we might expect each solution to contain two lines of symmetry, forming a 
rhombus. We show here that this is not necessarily true. When A' > and in G (m*, —1/2), 
where m* ~ —0.5951, there exists a family of convex kite central configurations that are 



not rhombi. There does exist a family of rhombi when A' > 0, as discussed in Section 7.4 
In fact, this family of rhombi undergoes a pitchfork bifurcation at m = m* that results in 
the convex kites. 
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Proposition 7.5. Let T 1 — T 2 — 1 and T 3 = T 4 = m < 0. Any strictly planar solution to 
the Albouy-Chenciner equations with A' = 1 must possess a line of symmetry. 



Proof. The type of computations required are very similar to those used in Sections 6.3 



and 7.2| Consider the Albouy-Chenciner and the unsymmetrized Albouy-Chenciner equa- 



tions together with the Cayley-Menger determinant and Dziobek equations with the nor- 
malization A' = 1. Take a Grobner basis of the ideal generated by such equations and 
saturate with respect to some of the variables to eliminate possible solutions where one 
of the mutual distances has zero length. Then, saturating with respect to the differences 
(S13-S24), O23-S14), O13-S14), O23-S24), O13-S23) and (S14-S24), yields a system of 12 
polynomial equations. Analyzing this system by the method of triangular decomposition 
of semialgebraic systems (with the inequalities Sjj > 0) yields an empty triangular system. 
Therefore, due to the saturation and the Symmetry Lemma, there are no asymmetric 
solutions in this case. □ 

To find the solutions with A' > 0, we impose the symmetry on the configuration and 
use equations (|7|) since the areas A4 are easily computable. Suppose that vortices 1 and 2 
are on the axis of symmetry and set r^ = r^ and r24 = ^3. We impose the scaling r^ = 1 
rather than choose A' = 1. Introduce two new variables x and y measuring the distance 
between the intersection of the two diagonals and vortices 1 and 2, respectively. Thus, we 
have the simple geometric equations 

x 2 + i = r 2 3 , y 2 + - = rf 3 and x + y = r 12 . (31) 



Equations (31) taken together with r 34 = 1 imply that the Cayley-Menger determinant 
ecu vanishes. 

Under this setup, convex configurations correspond to x > and y > 0, while concave 
configurations have xy < 0. The oriented areas of the four triangles are given by A\ = y/2, 
A2 = x/2 and A3 = A^ = —rw/4:, where the signs are taken without loss of generality. 
Note that one of the Dziobek equations in (|8| is automatically satisfied. The other equation 
yields an expression for A' given by 

22 2 

yl '13' 23 '12 



?12(?13 + r 2 3 ) - r l3 r 2 3 ( r l2 + 1) 

Equations ([7]), the formula for A' and equations ( |3~T| ) yield a polynomial system in the vari- 
ables r 12 , r 13 , r 2 3, x, y, X' and m. We then saturate this system with respect to r 13 , r 2 3, r 13 — 
r23 and r 13 + r 2 3 to eliminate the rhombus solutions and insure the mutual distance vari- 
ables are nonzero. Denote the resulting polynomial system as Pkite- Analyzing Grobner 
bases of the ideal generated by Pkite with different elimination orderings yields the following 
theorem. Recall that m* m —0.5951 is the only real root of the cubic 9m 3 + 3m 2 + 7m + 5. 

Theorem 7.6. For m* < m < —1/2, there exists four convex kite configurations with 
vortices 1 and 2 on the axis of symmetry. These solutions have A' > and are not 
rhombi. There are no other strictly kite solutions (with vortices 1 and 2 on the axis of 
symmetry) with A' > and m < 0. 
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Proof. Let Gkite be the Grobner basis of the ideal generated by P kite with respect to the 
lex order where r 23 > r 13 > r 12 > A' > y > x > m. The first element in Gkite is an even 
8th-degree polynomial in x with coefficients in m. Letting z = x 2 , this polynomial is given 
by 

C m (z) = 256m 2 (m + 2) (2m + l) 2 z 4 - 256m(9m 4 + 23m 3 + 17m 2 - m - 3) z 3 
+ (1728m 5 + 3136m 4 + 992m 3 - 384m 2 + 64m + 128) z 2 
+ (-432m 5 - 336m 4 + 48m 3 - 80m 2 + 16m + 64) z + (m + 2) 3 . 

Computing a Grobner basis with the same ordering of variables except that x > y > m 
yields the same polynomial with z = y 2 . The discriminant of ( m is a positive multiple of 

m 4 (m + l) 3 (2m + l)(9m 3 + 3m 2 + 7m + 5)(9m 2 + 4m - l) 2 (m - l) 2 (q k (m)) 2 

where gfe(m) = 2m 5 + 8m 3 + 14m 2 + 4m — 1. Note that the discriminant vanishes at m = m*. 
It also vanishes at m = —1/2, where Cm(z) becomes a cubic polynomial with roots 3/4 
and —1/4 (multiplicity 2). It is straight-forward to check that the discriminant is strictly 
negative for m* < m < — 1/2. Consequently, ( m (z) has precisely two real roots for each 
m-value in this interval. 

To see that these roots are always positive, we first use resultants to compute the 
repeated root of Cm(z) when m = m*. This yields the value z* ~ 1.9566, the only root of 
the cubic c u (z) = 320z 3 — 656,2 2 + 60z — 3. Next consider the ideal in C[z,m] generated 
by the system {( m ,c u }, which has (z*,m*) in its variety. Computing a Grobner basis 
for this ideal which eliminates z, we obtain a polynomial in m that has no roots for 
m* < m < —1/2. Since ( m (z*) < when m = —1/2, it follows that ( m (z*) < for any 
m G (m*,— 1/2). Then, since the constant term and leading coefficient of ( m are both 
positive for m* < m < —1/2, the two real roots of ( m must be positive. 

Choose x to be the positive square root of one of the roots of ( m . The second el- 
ement in the Grobner basis Gkite is linear in y, and the coefficient of y is non-zero for 
m* < m < —1/2. Thus, by the Extension Theorem, we can extend a solution (x,m) of 
Cm(x) = to a unique partial solution (y,x,m), where y must be real. By considering 



the other polynomials in Gkite, three of which are given by equations (31), the Extension 
Theorem repeatedly applies to extend (y, x, m) to a unique full solution of system P kitc 
where x, r 12 , r 13 , r 2 3 are all positive and r 13 ^ r 23 . As expected, many of the elements 
in Gkite featuring the saturation variable for r 13 — r 23 have a leading coefficient that is a 
multiple of the key cubic 9m 3 + 3m 2 + 7m + 5. 

To see that A' is positive, we compute a Grobner basis for Pkite which eliminates all 
variables except A' and m. The first element in this Grobner basis, a quadratic polynomial 
in A' with coefficients in m, is 

2m 2 (m + 1) A /2 + (4m - l)(m + l) 2 A' + m(m + 2)(2m + 1). (32) 

It is straight-forward to show that the roots of this quadratic are real and positive for 
m* < m < —1/2. Consequently, A' > and the configuration corresponding to the 
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solution guaranteed by the Extension Theorem is convex. Therefore, y > is also assured, 
and by symmetry, y must be the positive square root of the remaining real root of £ m . 
Since we have saturated to eliminate the rhombus solutions, we must have x 7^ y. Thus, 
there are two solutions: one where x > \fz* > y > 0, and one where y > \/z* > x > 0. 
The full count of four solutions comes from reflecting each of these solutions about the 
axis of symmetry (or simply interchanging vortices 3 and 4). This proves the first part of 
the theorem. 

To see that there are no other solutions (other than a rhombus) with A' > 0, first 
consider the case m G (—1/2,0). According to equation ( [32] ) , the two possible choices for 
A' are real, with opposite sign. Denote the larger, positive root as A' + and suppose there 



was a real solution to system P kite with A' = A' + > 0. Since the quadratic (32) is strictly 
negative when A' = 3 and m G (—1/2, 0), we can assume that A' + > 3. Next, the top two 
equations in (J7]) imply that 

im 2 xy{l + X') = 1 + AV 2 2 . (33) 
If our real solution has A' > 0, then it is convex and xy > 0. Since A' + > 3, we have that 

1 + X' + (x 2 + y 2 ) + xy(X' + - 1) > 0. 
Taking m G (—1/2, 0), this in turn implies that 

Am 2 xy(l + A') < xy(l + A') < 1 + X'Jx + y) 2 , 



which violates equation (33). Therefore, any real solution for m G (—1/2,0) must have 
A' < 0. 

Recall that for m = —1/2, ( m reduces to a cubic polynomial with only one positive 
root at ( = 3/4. Therefore, the only possible solution to system P^te at m = —1/2 has 
x = y = \/3/2, a rhombus solution. At m = m*, ( m has only one positive root at z* 
(multiplicity two), which also implies that x = y. For the case — 1 < m < m*, the 



discriminant of ( m is positive; however, using Lemma |3.5| it is straight-forward to check 
that all the roots of ( m are complex. Hence there are no real solutions to system P^te 
when — 1 < m < m*. 

For m < —1, the discriminant of ( m is negative so there are two real roots. To 
determine their sign, we note that C m (— 1/4) = 128m 3 (m+l)(2m+l) is strictly negative for 
—2 < m < — 1. Since both the leading coefficient and the constant term of ( m are positive 
when —2 < m < — 1, it follows that the real roots of ( m are negative in this case, and 
there are no real solutions to system Pkite- For m < —2, the signs of the leading coefficient 
and constant term both flip to become negative. Since Cm(3/4) = — 128(m — l)(2m + l) 2 
is strictly positive when m < —2, it follows that ( m has two positive real roots. However, 



examining the roots of the quadratic in equation (32 ), we see that A' < whenever m < — 2. 
Finally, when m = —2, ( m reduces to a cubic with only one real root at zero. Thus, we 
have shown that for m < and m (jL (m*, —1/2), there are either no real solutions or only 



solutions with A' < 0. This completes the proof of Theorem 7.6 □ 
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Remark. 1. Since Theorem 7.6 applies for all m < 0, we do not need to consider the 
case with vortices 3 and 4 on the axis of symmetry. If there existed a solution for some 
m G (—1,0) having A' > and vortices 3 and 4 on the axis of symmetry, we could 
rescale the vortex strengths by 1/m and relabel the vortices to obtain a solution with 
m < — 1, A' > 0, and vortices 1 and 2 on the axis of symmetry. But this contradicts 
the last statement of the theorem. 

2. For completeness, we note that the cases having real, positive solutions with A' < 



agree with the results in Section 7.2, For example, when —1/2 < m < 0, we find 



concave kite configurations with vortices 1 and 2 on the axis of symmetry, as predicted 



by Theorem |7.2| For m < —2, we find convex kite solutions with vortices 1 and 2 
on the axis of symmetry. Rescaling the vortex strengths by 1/m and relabeling the 
vortices gives a family of convex kites with vortices 3 and 4 on the axis of symmetry 



and — 1/2 < m < 0. This also concurs with Theorem 7.2 



7.4 The rhombus solutions 

In this final section we study convex kite configurations that have two lines of symmetry, 
namely rhombus configuations. Here we have four congruent exterior sides, 7*13 = ru = 
r 23 — r 24 with the diagonals satisfying the relation r\ 2 + r| 4 = 4r 2 3 . The areas satisfy 
A\ = A 2 = —A 3 = —A4. The first and last equations in (|9j) are clearly satisfied. The 
middle four equations in ^ are equivalent. They, along with the value of A' determined 
by (18]), yield a formula for m given by 



- = ^> (34) 

where x = r^/r 12 is the ratio between the diagonals of the rhombus. Note that m > if 
and only if 1/a/3 < x < \/3. As usual, we restrict to the case — 1 < m < 1. 

Unlike the isosceles trapezoid case, here we have solutions for m < 0. In fact, there are 
two geometrically distinct families of rhombi when m < 0. This can be seen by inverting 



equation (34), which yields 

- (3 - 3m ± v/(3 - 3m) 2 + 4m ) , (35) 



x 2 

2 



where + is taken if m > and + or — may be chosen if m < 0. Using equation (El), we 
compute that 

A = - 4K + 4m + l) 

rj 2 (2 + 3m — 3m 2 ± my 9m 2 — 14m + 9 ) 



with the same sign choices as for x . Note that the numerator of (36) vanishes at m 



-2 + \/3 ~ —0.2679. The denominator also vanishes at this special value, but only when 



+ is chosen. For the rhombus solution when + is taken in both (35) and (36), the value 



of the angular velocity A is always positive (thus A' < 0) and monotonically increasing in 



m for m G (—1, 1]. In contrast, for the rhombus solution when — is taken in both (35) 
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Figure 4: A plot of x versus m for the two rhombus families, where x 
between the diagonals of the rhombus. 



^34/^12 is the ratio 



and Q36J) (m < only), we have A > (thus A' < 0) only for -2 + a/3 < m < 0. At 
the special value m = —2 + \/3, the rhombus relative equilibrium actually becomes an 



equilibrium (see Section 4.1). Asm decreases through —2 + \/3, the direction of rotation 
flips, A becomes negative and A' becomes positive. We summarize our conclusions in the 
following theorem. 

Theorem 7.7. There exists two one-parameter families of rhombi relative equilibria with 
vortex strengths I\ = T 2 = 1 and T 3 = T 4 = m. The vortices 1 and 2 lie on opposite sides 
of each other, as do vortices 3 and 4- Let (3 = 3 — 3m. The mutual distances are given by 



fe^j = I (/?± v//32 + 4m) and (^\ = ^ (/3 + 2 ± a//3 2 + 4m 



(37) 



describing two distinct solutions. Taking + in (31) yields a solution for m e (—1, 1] that 
always has A > 0. Taking — in (31) yields a solution for m G (—1,0) that has A > for 
m G (-2 + a/3,0), but A < for m G (-1, -2 + y/3). At m = -2 + y/3, the - solution 
becomes an equilibrium. The case m = 1 reduces to the square. Form > 0, the larger pair 
of vortices lie on the shorter diagonal. 

Proof. The formula for r^/r 12 comes from 1 + x 2 = ^{riz/r^) 2 . For the case m > 0, it is 
easy to show from equation (34) that m < 1 if and only if 1 < x < y/3. This verifies the 
last statement in the theorem. □ 



Remark. 1. One way to visualize the two rhombus solutions is captured in Figure |4j 
where a plot of x = r^/r^ versus m is shown. Beginning with the square at m = 1, 
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Figure 5: The two distinct rhombi relative equilibria when m = 
in opposite directions. 



—0.3. The solutions rotate 



as m decreases, the ratio of the diagonals of the rhombus increases from 1 to \/3, with 
the larger pair of vortices lying on the shorter diagonal. At m = a bifurcation occurs, 
and a new family is born emerging out of a binary collision between vortices 3 and 4. 
Unlike the previous family, which continues past the bifurcation, this family has the 
larger pair of vortices (in absolute value) on the longest diagonal (see Figure [5]). As 
m — > — 1 + , x — > v2 — 1 for the new family while x — > v2 + 1 for the previous family. 



2. The rhombus solution when — is chosen in equation (37) undergoes a pitchfork bifur- 
cation at the special parameter value m* pa —0.5951. As m increases through m*, the 
rhombus solution bifurcates into two convex kite solutions with the positive strength 
vortices on the axis of symmetry. The two kites are distinguished by whether ri3 > 
or ri3 < T23. Since the rhombus solution continues to exist past the bifurcation, we 
have a pitchfork bifurcation. The Hessian matrix D 2 (H + XI) evaluated at the — 
rhombus solution at m = m* has a null space of dimension 1 (excluding the "trivial" 
eigenvector in the direction of rotation) and contains an eigenvector corresponding to 



a perturbation in the direction of the convex kite solution found in Section 7.3 



3. The fact that there are two geometrically distinct rhombus solutions for m < (plus 
two convex kites when m* < m < —1/2) indicates that there is not a unique (up to 
symmetry) convex central configuration for these four choices of vorticities. This con- 
trasts with the Newtonian four-body problem where it is thought (although unproven) 
that there is a unique convex central configuration for any ordering of four positive 
masses. 
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